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TESSERAL MATRICES 
By T. SMITH (Teddington) 
[Received 15 June 1931] 


1. WHEN it is desired to exhibit the number 7 of rows and the 
number p of columns of a matrix u the notation _w? will be employed. 
The orders of a product are then given by the external subscript and 


superscript: thus 
I w= Uv = wv? = wv’. 


Suppose that p, r, s are definite non-singular square matrices of orders 
7, p,o and that P, R, S are their inverses. From wu, v, w we form 
new matrices U, V, W, 
U = puR, V=rv8, W=pwS, ¢ 

by prefixing the matrix corresponding to the order of the subscript 
and postfixing the inverse matrix corresponding to the order of the 
superscript. Then 

W = puvS = puRrvS = UV; W’ =7'U", 
where accents denote conjugate matrices. It is evident that, if U’ 
and V’ are the inverses of uw and v, then W’ will be the inverse of 
w provided the sequence 7, p, o is monotonic. 


2. Now suppose that the relation between u and U is symmetrical, 


so that u = pUR = puR*. 

This equation will hold for a large class of matrices if p? and R? are 
scalar, or more precisely if p = «P, r = «R, where « is a definite scalar. 
No loss of generality results from assuming |«| = 1. If « = 1, our con- 
dition leads to a group of which orthogonal matrices are the simplest 
examples. If « = —1, the simplest class is obtained by taking 


=i 
Pa SS 


where 0 and 1 denote zero and unit matrices of the same order. The 
orders of u are then necessarily even, and it is convenient to replace 
7, p, o by 2m, 2p, 2o. A typical matrix wu of this class is related to 
| its reciprocal U and its inverse U’ as shown by the equations 


, 


a , fl) 
=(0 > 8=(a 6 =(_o oe) @ 
, l U’ | (2.1) 


where a, b, c, d are each of the form ,q?. I shall refer to members of 


‘this class as tesseral matrices. 
8695.2 R 
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Tesseral matrices are necessarily non-singular, and the way in 
which they have been derived shows that the product of tesseral 
matrices is a tesseral matrix if the orders form a monotonic sequence. 
3. In what follows we assume that at least one of the components 
a, b, c, d is non-singular. This is not a general property of tesseral 
matrices. For example, 


] 


is a tesseral matrix having a, b, c, d all singular. 


The four components a, 6, c, d are not independent, and the con- 


nexion can be shown explicitly. Consider the case when a is non- 
singular and 7 <p. The appropriate condition wU’ = 1 implies 
bc’ —da’ = 1, cb’—ad’ = 1, ac’ = ca", bd’ = db’, 
three independent relations, for the second is the conjugate of the 
first. The third condition can always be satisfied by 
c= ey, ” == fa’, 
where f is symmetrical of order p, and then the first condition if 
d= bf—A, 
where A is any* solution of Aa’ = 1. Substitution in bd’ = db’ gives 
) A’ = Ab’, and the general solution is 
bA 1] g 
b= AE, b’= EA’, 
where E is symmetrical. 
The resulting expressions for u and U’ are 

iw) 

AK AEf— * v’ fa’ A’ Cg 

a af -a' EA’ }’ 
from which it is clear that the conditions wU’ = 1, U'u=1 reduce 
to aA’ = 1, A’a= 1 respectively. 

We shall be concerned only with the slightly less general solution 
derived from the substitution 
b= ea, b’ =a'e. e= AEA’. 

For convenience we shall call tesseral matrices of this type t-matrices. 
The four forms for 7-matrices corresponding to non-singular values 


* A is only determinable uniquely from @ when 7 = p. 
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of a, b,c, d are 


eal PI” 2 


) Say 
om is bh r= (Were es | 24) 
= 


fa’ eae ie | (3.1) 


gb = =gbh+-B)’ —b'g b’ 
Se a Sel *) c’ 


) 
cj c ae jcei+t ae 5 


a 4 
=e uw Ue bites m4 (3-4) 


If += 4,77", then e, g, i, k are symmetrical of order z and f, h, 7, 1 
symmetrical of order p. The condition 77’ =—1 depends upon 
aA’ =bB' =cC’ =dD' = 1 in these respective cases, and similarly 
T's =1 upon A’a= Bb = C’c= D'd=1. When any of the com- 
ponents a, b, c, d is singular the corresponding 7 form does not exist. 
No confusion will arise, if we also regard w in (2.1) as a 7-matrix, 
provided it can be put into one of the forms (3.1), (3.2), (3.3), (3.4). 


4. We now proceed to extend the multiplication theorem by show- 
ing that the product of 7-matrices is a 7-matrix subject only to the 
obvious conditions concerning the rank of the product components. 
It will be sufficient to take the A form of 7 as representative of the 
A and D forms, and the B form as representative of the B and C 
forms. Also we need only take the cases 7 <p <o and <p Do, 
= <o to represent the four cases that may arise. By direct multi- 
plication we have 


eh ee oy “et, oe 


a, a,f; ay asf. a af 
identically, where a = a,(f,+¢,)de, (4.1) 
(e,—e)a = Ajay, a(fo—f) = a,Ag. (4.2), (4.3) 


The components of the product matrix are linear in the magnitudes 
represented by the compound symbols ea and af. We have to show 
that it is legitimate to represent them in this way as products, that 
e and f may then always be symmetrical, and that (e,—e)a(f,—f) is 
reciprocal to a. The notation is clearly justified, as (4.2) and (4.3) 
can be solved for e and f. 
If <p <a, the equations immediately yield 
R2 
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aAsa’ = a,(f,;+e.)a,, a symmetrical matrix, = a,A,0’, 
(e,—e)a Aja, = 1, a(f.—f)a’ = a,A,q’. 
These conditions require e to be a symmetrical matrix of order 7, and 
always permit f to be symmetrical of order o. Moreover, 
(e,—e)a(f.—f)a’ = 1, 
so that we are justified in writing 
(e,—e)a(f.—f) = A, 
which completes the proof for this case. 


When z <p ><, put 6 =a,A’a,. (4.4) 


From (4.1) we have 
6(f,+e,)0 = a,A’aA’a, = 8, (4.41) 
when either aA’= 1 or A’a=1. This equation always admits of a 
symmetrical value of 6. When @ is known, A’ is found from 
A = ADS, 
and 
(€;—e)a(fo—f) = (€,—e)aA’a(f,—f) = Aya,A’a,A, = AOA, 
If now 7 <a, 
aA’=1 and e,—e=A,a,A’=A,0A);, (4.43) 
so that e is necessarily symmetrical, and 
a(f,—f)a’ = aAj0A,a’, (4.44) 
showing that a symmetrical form for f exists. This condition is clearly 
met if we take the solution 
f.—f = A,0A, = A’a,Ayz. (4.45) 
When z = p = 9, (4.43) and (4.45) represent the only solution which 
exists, and it is clear that they are always acceptable solutions of 
(4.2) and (4.3). We remove all ambiguity by invariably adopting 
these expressions for e and f. 
Turning now to the form of 7 involving B, we have 
b, b,h, \ b, bohy 


fx, Gxbyhyt+ By \ Gob. Jaboha+ B, 


b bh 
(1 gbh+- B,B,—(g—g,)b(h—h,) ): 
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where 


b = b,(1+hyge)bo, 
(9—9)b = By gabe, 
b(h—hg) = b,h, B,. 


Taking first the case 7 < p <a, (4.5) gives 


b Bihyb, = b,(h, +h, goh,)b;, a symmetrical matrix, 
= bh, B,b’ = b(h—h,)b’, 
so that a symmetrical value of h is possible. Also 
(g—g,)b(h—h.)b’ = By, goh,b, = B,(goh,+1)b,—1 = B,B,b’—1, 
which justifies us in writing 
B= B,B,—(g—g,)b(h—h,), Bb’ = 1. 
From this equation 
Bb; Job2 = [ By—(g—gy)byhy gabe 
= (9—91)(6—b,hy gob) = (9—91)byb, 
or Bbi,qB; = 9—g, — By gaboB, 
so that g is necessarily symmetrical. 
Coming now to the case 7 < p >a, put d= 6, B’b,. Then 
dh = b, B'bB'b, = d(1+h, ge), (4.8) 
showing that values of ¢ can be found which give symmetrical values 
of dh, and g,¢. It is not to be assumed that both are symmetrical 
for the same value of ¢. Take g,¢ = ¢’g2, and multiply the equation 
for ¢ behind by h,¢’. Then 
7h, bd’ = th, d’—h, .go¢-h, ¢’, a symmetrical matrix, 
= gh, ¢”’, 
from which we can justify the equality 
gh, =h,¢’. 
Also, if™<o, — g_g, — By gyb,B’ = B,g, $B; 
and b(h—h,)b’ = bB, dh, Bp’, 
so that g is necessarily symmetrical, and a symmetrical value can be 
found for h. If we now assume that B satisfies the equation 
B, = Bby1+-geh,), 
we note that the conditions of order and rank are satisfied. Multi- 
plication behind by 6; gives 


Bb, — Bb’, 
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by (4.5), multiplication in front by b; and behind by ¢’ gives 
d' = $'(1+goh;)¢’ 

in accordance with (4.8), and B, 9,b,B’, i.e. g—g,, is obviously sym- 
metrical. We are therefore justified in adopting this equation as a 
definition of B. Now multiply the equation behind by B,. Noting 
that Bb, = B, 4’, we have 

B,B,— B= B,d'g.h,B, = B,92.¢h,B, = (g—g,)b(h—h,), 
which completes the proof of the theorem. 

As with the A form, we may always render the solution of (4.6) 
and (4.7) unambiguous by confining our attention to the particular 
solution g—9, = B, gobB’, (4.91) 

h—h, = B'b,h, By. (4.92) 

5. When more than one form for a particular r-matrix exists, rela- 
tions must exist between their elements. Let us suppose 7 < p. 
Comparing the first columns of (3.1) and (3.2), we have 

a=4@), b= ea, a= gea, 
or ge=—1=@, 
on multiplying behind by A’. 
Again, from the first columns of (3.1) and (3.3) 
ea = 1a+C, 
or e=1+CA’. 
Finally, from the second columns of (3.1), (3.2), and (3.4), 
A = eaf—d = e(gbh+ B)—d = eB. 


These relations are typical of a complete set involving e, g, 7, k, viz. 
A =eB, B=gA, C=wW, D = kC, 
e—i=CA'=AC’, k—g= BD'= DB’, 
eg = ge=1, k= ie = 1, 


which also involve 
ek—AD’ = 1, git BC’ = 1. 
We next obtain equations connecting e, g, i, k with f, h, 7, 1. From 
(3.1) and (3.4) 
Al = (ec—d)l = e(D+a)—dl = eD (5.51) 
and fD' =f(b’'k—a’) = (A'+d')k—c' = A’k. (5.52) 
Similarly gC = Bj, BY =hC". (5.53) 


Finally, considering the relations involving f, h, j, 1, we obtain equa- 
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tions corresponding to (5.1) 
Cf = (b—1a) f= A+d—ic=A 
directly, or alternatively by (5.1), (5.3), and (5.51), (5.52), (5.53), 
A= G, B= Dh, C= Aj, D= Bil. (5.6) 
When 7 >p we may also establish relations corresponding to (5.2), 
(5.3), and (5.4), but when 7 <p we can only infer lower-order rela- 
tions such as a(fj—1) = 0, owing to the failure of a and af to specify 
f uniquely. When the ambiguities in the symmetrical factors of order 
p are avoided as in the last section, the limiting factors such as a can 
be removed and one form of the 7-matrix can be obtained from 
another. All the remaining equations can be derived from a suitable 
set of three. For example, noting that 
A’ B = (c’e—d’)(ec—gd) = (c’—d’g)(ec—d) = B’A 
by (5.3), and similarly C’D = D’C, we see that the assumptions 
f—h=A'B= BA, l-j = C’D=D'C (5.7) 
are consistent with the symmetry of f, h, j, 1. The first of these 
assumptions involves the equivalence of (3.1) and (3.2), and the 
second that of (3.3) and (3.4). To complete the connexion of the 
whole series, we assume fj = 1. From (5.7) 
fG+C'D) = (h+A’'B)l 
or fji+A’'D=hl+A’'D 
by (5.6), so that we have both 
5-7 — 1, hi=th=1. (5.8) 
By (5.6), (5.7), (5.8) we easily derive 
fli—A’D=1, hj+ B’C=1. (5.9) 
6. The elements of the standard z-matrices (3.1), (3.2), (3.3), (3.4) 
which are absent from the u-matrix (2.1) may be arranged to form four 


symmetrical square matrices of order 7+-p. The definitions are 


“=(_4 5) 8=( 4)] 


D C s 2 
= $= ; 
? (o _;) ' ES L, 


(6.1) 


It will be noticed that the symmetry of a and 6 is more complete 
than that of B and y. This may be ascribed to the displacements 
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which have to be imparted to the elements of the second and third 

forms of 7 to bring them into the positions they assume in f and y. 
By (5.4), (5.51), (5.52), (5.53), (5.9) we see that 


m= Sa, By = 1= yf. (6.2) 


7. We now consider the bearing of the relations that have been 
found on contact transformations. In the first place we shall deal 
with the significance of the 7 multiplication theorem, taking as an 
illustration the A form associated with the «-matrix, and afterwards 
extend the result to the other forms by changes of variables. Let 
x, y, z be matrices of variables having the same number of rows and 
7, p, « columns respectively. As before, let the orders of 7, be 27, 
2p, of 7, be 2p, 20, and of 7 be 27, 2c. 

Now the conditions that identically 

x’ y’ x 
(a y)ox ( ,) + (Y Z)os A (a z)a( ,) 
y’, z z 
= (wA,—p+zA)5)0(Aja2’—p’+ 
y = pe, O= ff, 
6 = O(f,+e.)8, A= A,0A,, 
e= ASA. fo—f = Aj0Az. 


€} 


The last four equations correspond exactly to (4.41), (4.42), (4.43), 
(4.45), where we also assumed that (7.3) was satisfied. Thus, when 
7 <p ><a, the left side can be factorized as shown on the right, pro- 
vided the magnitudes of the elements in « are those which would be 
derived from the r-equation 7 = 7,75. 

If we substitute from (7.2) and (7.31) in (7.1), we may write the 
equation in the form 


, 


x y’ x 
(x y)a(",)+y 2)a2(%,)—(@ aa(* | 
y z z 
(wAa,—y+z2A'a,)(f,+6)(a,A’x’ —y’+a)A2’), (7.4) 
where the value of @ given in (4.4) has been inserted. 
When 7 < p<o we consider the conditions that 


(@ Y)oy (*") +(Y Z)ote (",) (x n(®,) 


- [w(e,—e)a—qa+-z|(f.—f)| a’ (e,—e)a’ —a’q' +2’ | 
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should be an identity. As before, matrix multiplication of the 7’s is 
indicated, together with y = ga,. Since 

(€,—e)a, = (€,—e)a(f.—f)a, = Aay, 


the identity may be expressed in the form 


(x ya (J+ z)ag(¥, )— (x )n(*,) 
= (wAa,—y+2zA'a,)Aja(f.—f)a’A,(a,A'x'—y'+a\Az’). (7.5) 
Now since the extreme factors in the products on the right sides of 
(7.4) and (7.5) are the conjugates of one another, the left sides and 
their differentials with respect to y vanish when 
y = rAa,+zA'a,. (7.6) 


Consider, then, the problem of evaluating the sum 
2 2 ‘y’ 
(x y)a J+ z)a( ,). 
"\y 2 


where a, and a, are symmetrical matrices, and y is to receive the 
value which makes the sum stationary for any small variations in y. 
In general «, is a function of x, y, and a, a function of y, z. Obviously 
the result can be expressed in the form 


where « is a symmetrical matrix in general depending on x and z. 
The equations (7.4) and (7.5) show that « is determinable from the 
matrix T= 7,72, provided the value of y given by (7.6) is inserted in 
the product. This equation for y is convenient for solution by suc- 
cessive approximation, and is very suitable for use with the matrices 


since the effective functions a, and a, are known. This form of solu- 


tion has obvious advantages over that derived by a direct attack 
upon the problem. 

8. Let x, y, €, y be matrices which at the moment are only restricted 
to having the same number of rows and 7, p, 7, p columns respectively. 
For convenience write also 
k=(ry), A=(€ -y), w=(@—n), v= —7) (8.1) 
and A= kak’, B = Ap’, C= py’, D = viv’. (8.2) 
By (6.2) we have identically 

(ea+v)(«’+8v’) = A+D+xn0'’+rx’ (8.3) 
and (AB -+1)(0 + yp") = B+O+Aw'+ pr’. (8.4) 
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Also, by various relations given in § 5 
(xe—yA'+€&)(x’— By’ +g’) aa | B+ wé’+ £2’ (8.5) 
(ai+€—7C’)(x’ +ké’— Dn’) = D+ C+22’ a (8.6) 
(xA —yf+)(C’x’—y’ +jn') = A—C—yn'—ny’, (8.7) 
(yht+EB—n)(y'+D'e—ly’) = D—-B —yy' —ny’. (8.8) 
Now in all these six identities the two factors on the left are obtain- 
able from one another by transposition and multiplication by a 
suitable factor. It follows that by making these factors vanish we 
obtain contact transformations involving the expressions on the right 
of these equations. It is easily verified that x, y, €, » are used con- 
sistently in these equations. Thus all the twelve left-side factors 
vanish if cxa+v = 0 in virtue of the relations of § 5. We are therefore 
able to construct the remaining functions of the group A, B, C, D 
from a knowledge of any one of them, provided more than one exists. 
For example, given A as a function of x and y, we can determine 


D a function of £ and 7» defined by 


, 


, 


A+D-+ 10’ +k’ i dD 0 


dk 


by substituting from «+v5= 0 in 6. Transformations of this type 


exemplify the principle of duality. 

The slightly imperfect symmetry of 8 and y appears again in the 
transposition of the Greek and roman letters in passing from A to p. 
The negative signs in the second members of A, u, v depend upon the 
definitions adopted for 8, y, 5. The signs in these expressions have 
been chosen for convenience and to enable the standard relations to 
include the simple equations (6.2). 

To the foregoing relations we add an equation of a different type 
which is involved in «+6 = 0, viz. 

(y —) = (x —§)r (8.9) 
and its inverse (y —n)T” = (x —6&). (8.91) 
This equation is suitable for use when 7 is given initially as a function 
of x and €. It shows that the problems which give rise to these 
methods can be solved in a different way. It should be remembered 
that the relations which subsist between the elements of tesseral 
matrices only involve their magnitudes, and a particular set of magni- 
tudes which is here always represented by a single character may in 
different equations be different functions of different variables. 
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9. Expressions which involve both the small and capital roman 
letters, such as (7.6), are formally indeterminate when the com- 
ponents are not square. The difficulty is almost always only apparent. 
In this example we may adopt 2A and zA’ as the effective variables. 
It will almost invariably happen, however, that the final product of 
the 7’s is a square matrix, and there is no ambiguity in the expres- 
sions ultimately obtained. The initial conditions also usually give 
rise to square matrices. Components with unequal orders are fre- 
quently the result of simplification in the early stages of a problem. 
Some elementary results of importance may be noted. 
(i) If b =c=1 and either a = 0 or d = 0, the matrix is a r-matrix. 


1 a, Ce, 
the same side with their inverses to yield scalars, so is 


(ii) If (" _ and % “:) are tesseral matrices multiplying on 
\% © 


; Cy 
and if the first two are r-matrices, so is the third. The converse of 


this proposition is also true. 


(iii) If (, “) is a tesseral matrix, so is the matrix obtained by 
a 


deleting any corresponding long rows* from a, b, c, d. Since the 
distinguishing property of a 7-matrix is that the long rows of at least 
one of the four components a, b, c, d are independent, we see that the 
contracted matrix is a 7-matrix, if the original matrix is a r-matrix, 
but the converse of this proposition is not true. 

(iv) The deletion of any long rows from a, b, c, d will be accom- 
panied by the deletion of the corresponding rows from A, B, C, D. 

(v) Only two components of a tesseral matrix may have corre- 
sponding long zero rows, and these must be a diagonal pair. 

(vi) Corresponding short zero rows may be added.to or deleted 
from the four constituents a, b, c, d without affecting the character 
of a tesseral or r-matrix. 

(vii) In general, short rows may not be removed from tesseral 
matrices. 

* i.e. either rows or columns, the choice depending on which has the greater 
number of elements. The other set of rows are the short rows. 








A PROPERTY OF INTEGRAL FUNCTIONS OF 
FINITE ORDER 
By J. M. WHITTAKER (Cambridge) 
[Received 21 May 1931] 
1. A WELL-KNOWN theorem due to Littlewood, Wiman, and Valiron* 
states that an integral function of order less than a half is ‘flat’, i.e. 
log m(r) > a positive constant x log M(r), 


on a sequence of circles of indefinitely increasing radius. In a recent 
papery is was shown that these circles can be broadened out into 
annuli of flatness; in fact that 


log m,(r) > a positive constant x log M,(r) 


for a sequence of annuli r< |z| <r+r?’, (co<1—p), m,(r), M,(r) 


denoting the bounds of |f(z)| in the annulus. 

Functions of higher order need not possess annuli of flatness, but 
it was shown in I that in the case of a function of genus 0 or 1{ there 
is a sequence of circles of fixed radius within which the function is 
flat. Thus, 


THEOREM A. Let f(z) be an integral function of genus 0 or 1 and let 
d be a given positive number. Then there is a positive constant h and 
a sequence C,, Co,... such that 


A { |>a 


Ss 


and log |f(z)| > hlog M(A,) 


8 


in the circle |z—C,| <d. 

* Littlewood, Proc. London Math. Soc., 6 (1908), 189-204; Wiman, Math. 
Ann., 76 (1915), 197-211; Valiron, Ann. de Toulouse, 5 (1913), 117-257. Proofs 
have also been given by Polya, J. London Math. Soc., 1 (1926), 78-86; and 
Besicovitch, Math. Ann., 97 (1927), 677-95. 

+ Proc. Edinburgh Math. Soc., 2 (1930), 111-28. This will be referred to 
as I. Miss Cartwright has pointed out to me that Theorem 2 of I, on functions 
with zeros at the lattice points, can be deduced from a result of F. and R. 
Nevanlinna, ‘Uber die Eigenschaften analytischer Funktionen in der Umge- 
bung einer singuliren Stelle oder Linie’, Acta Soc. Scient. Fennicae, 50 (1922), 
p. 42. eo 

P , ( log M(x) 
i.e. such that J 3 da 
1 


converges. See Valiron, Integral Functions, Toulouse (1923), p. 185. 
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In the present paper a result of this type applicable to functions 
of any finite order is given. 
THEeorEeM ©. Let f(z) be an integral function of order p, and let o be 
a given number <1—4p. Then there is a positive constant h and a 
sequence L,, Cy,... such that (1) is true and (2) is true in the circle 
je—{,| <A¥. (3) 
This result is the best possible in the sense that o cannot in 
general be replaced by 1—$p in (3). 


Theorem C is suggested by writing z = w'? in Theorem A. Of 
course this procedure does not constitute a proof, but the demonstra- 
tion which follows can be looked on as a rigorous formulation of the 
same principle. The details are rather long and have been omitted 
for the most part. 

Lemma 1. Let f(z) be of genus 0 or 1, and let X, € be given positive 
numbers. Divide the z-plane into squares of side X by drawing lines 
parallel to the axes. Then for almost every square 

llog f(2,) — log f(zz)|< er(logr)', (4) 
where z,, 2, are any points of the square and r is the distance of the 
centre of the square from the origin. 

This is essentially Lemma 6 of I. The conclusion of the lemma 


means that Vv 


vy 0, as No, 


where N, is the number of squares, among the N nearest the origin, 
for which (4) is false. 


Lemma 2. Let f(z) be an integral function of the form 
Cz x (a canonical product of genus p) x exp(b,2+-b,2?+-...+-b,2%), (5) 
having no zeros in the sector ~<0<B, (z= re'®). Then, if 
a+2y <6 < B—2y and |h|<rsin 2y, 


« 


1 qd. 
logs — log fie) plh|r -  Tn(Tn +7) +H 2, saat 


n= 


- 


h|\p+2 a | h|p+s 
 cosee” Hy > —; +——; cosec?* Ds —+ 
Pr- 1 (r,+1r)?* p+3 er 


n 


(6) 
where A, ,, depends only on y, p and rj, 1p,... denote the moduli of the 


ZETOS A, Ag... of f(z) 
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We make use of the inequalities 


a,+2| >(r,+r)sin y, a,+z| >rsin2y>h, 
which are readily proved with the aid of a diagram. Write 
9 
‘ oe. uP 
E(u, p) = (1- uexp| U+—s+..+-— 
so that as P 


log B| a P| — log E(—~.p) 
h Bie. a 1(—) z- _(—yp zP 


p 
a n Pan 


a 


-log(1 { 


a,+z 
and, on expanding the logarithm and making use of (7), we arrive 
at (6). © 

Lemma 3. Let S 


hyperbolae 


mn be a mesh formed by the intersections of the 

x*—y? = m, Icy = 8, x?—y*? = m+, 2zy=n+1, (8) 
where m, n are positive integers, lying entirely within the domain De , 
defined by the inequalities 

x<O0<}r—o, 4Rr<R, (x+iy=re, 0<a<}n, RD>1). (9) 
Then there are constants A, B, C depending only on x such that 

(i) AR~* < area of S,,,, < BR, 

(ii) Bas 

With the aid of a diagram this is readily seen to be true. The 
hyperbolae are the transforms of the straight lines of Lemma 1 (with 
A= 1) when the transformation z2 = w is made. Similar results hold 
for the transformation z” = w where k is any integer. The hyperbolae 


contains a circle of radius CR-. 


are replaced by the curves 
rkcoskO=m, résinkO=n, r*coskO=m-+1, = résink@=n-+1, 
and the domain D,,, by the domain 


=; te ob. 


7 
ie = 0b. IR<r<R, O<a< 
2 7 ( ‘ce 


2k 
The meshes contain circles of radius C R!-*. 
3. We next consider the ‘surface-density’ of log|f(z)|. 


For a domain D of area A(D) this may be defined to be 


o(f,D) “| logifte)| dedy. 


) 


A D ) | 
I 
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[It is shown in I (4.4) that for a circle of radius R this becomes 


R 
o( R) = log|C| + slog + fre="(1- 3) du, (10) 


where n(w) denotes the number of zeros of modulus not greater than 
wu, and f(z) is given by (5). In the case of an annulus }R<r<R 
the surface-density may be written 


+(R) = —_ (7 R8o( R)—}n R°o(4R)}, (11) 


and, since, by (10), o() is an increasing function of R, 

o(R) <7(R). (12) 

Thus o(R) can be replaced by 7(R) in Theorem 5 of I, which is then 

as follows: 
Lemma 4. If f(z) is an integral function of finite non-integral order, 
" R) 

lim 7( - 

r+» log M(R) 

If f(z) ts of finite integral order, this property holds for the function 
{ f(z) +d} for every value of d except, possibly, one value. 


> 0. 


4. We have now collected enough material to prove a preliminary 
theorem. 

Let f(z) be an integral function of genus p, order p. It is convenient 
to distinguish six cases as follows, 1 denoting a positive integer or 
zero: 

(i) p= 2, p= 2. (ii) AW2W<p<2%+1. (iii) p= 2, p= 2+1. 
(iv) p=2i+1, p=2t4+1. (v) 2+1< p< 2142. 
(vi) p= 21+1, p= 2242. 

Lemna 5. If f(z) falls under one of cases (i), (ii), (v), (vi), then (2) 

is true in a sequence of circles 
lze—f,| < da," 


[t will appear later that the lemma is true in cases (iii), (iv), but 
the proof is different. 
Consider first case (v) with i= 1, so that 3< p< 4. Let 
R a+tr 


: log |f(we!®) |w dudd 


127 R? 
tk a 


T,(R) wre 
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denote the surface-density in the domain 
tR<r<R, a<b6<a+ hn. 

Then for at least one of the values 4:7, #7, }:7,... of « 

: R 

lim 72) 

Rx log M(R) 

1 


There is no loss of generality in taking this value to be 47. We 


(13) 


can write f(z) in the form 
f(2) =fh@fr(2), 


where f(z) is a canonical product with zeros b,, by,... all of positive 


real part, and f,(z) is an integral function with zeros c,, ¢,,... all of 


zero or negative real part. Write 
Sil fil—2) = 9), 
so that g(w) is an integral function of genus 1 of the complex variable 
w= u-+iv, with zeros at the points w = bj, b3..... 
Now take positive numbers A, « and divide the w-plane into squares 
of side A. Shade the squares for which the inequality 
log g(w,)—log g(w.)| < et(log t)! 
is true; w,, w, are any points in a square whose centre is distant ¢ 
from the origin. 
Next divide the sector 0<@<}m of the z-plane into meshes 
formed by the intersections of the hyperbolae 
x?—y? = ma, 2ry=nrA, x?—y?=(m-+1)d, 2ry = (n+1)dr. 
and shade those meshes which correspond (by the transformation 
2? = w) to shaded squares in the w-plane. In a shaded mesh 
log g(z7) — log g(z3)| << 2er*(log r)* (r* = £), 
and, applying Lemma 2 to the functions f,(z), f;(—z) and remembering 
that |z,—z,| = O(r-), it is found that 
log fo(z,) — log f2(z2); = O(r*), 
a similar inequality holding for f,(—z). Now 
f(z) = 9) falz)/fi(—2), 
so that in a shaded mesh 
log f(z,) — log f(z,)| < 3er?(log r)! (r > 1). (14) 
Now it has been shown that there is a constant A, such that 
7,(R) > K, log M(R) (a = 777), (15) 


for a sequence f,, R,,... of values of R. 
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By Lemma 3 the meshes in the domain 
bR,<rcR, fet SOS fez 

are, roughly speaking, of the same size, and by Lemma 1 almost 
every one is shaded. The argument can now be completed as in the 
proof of Theorem A in I. Since, by (15), the average value of log|f(z)| 
is of the same order as log M(R,), it follows that there must be at 
least one point z, in a shaded mesh such that 

log f(z;) >Kzlog M(R,) >K,Re” — (p(B) >3), 


and so by (14) the inequality 


K, 
log f(z.) > log M(R,) 


must be satisfied for all points z, in the mesh. Finally, given d, it 
follows from (ii) of Lemma 3 that A can be so chosen that each mesh 
contains a circle of radius dR;'. 

This completes the proof of case (v). Cases (i), (ii), (vi) for i= 1 
can be proved in much the same way. 

The argument can be extended to prove the lemma for 7 = 2, 3, 4,.... 
The appropriate transformation is z'+! = w, and Lemma 3 is replaced 
by the extension of it involving the curves 

r'+1 cos(i+1)0 = m, ré+1 sin(t+1)0 = n. 
The auxiliary function g(z) is defined by 
gz") = f,le)fy(w2)...f,(w"2), 
where 1, w,...,w' are (i+1)th roots of unity and f,(z) is a canonical 
product having the same. zeros as f(z) in the sector 
7 7 
—.—- <9<._.. 
+1 +1 

5. We can now enunciate a more precise result* of the same 
character. 

THEorREM B. Let f(z) be an integral function of genus p. Then (2) 
is true in a sequence of circles 

pL] < dag 

As before, d is any given positive number. 

As regards cases (v), (vi) this has been established in Lemma 5. 

* This result has been made considerably sharper in accordance with a 


suggestion of Miss Cartwright. 


3695.2 
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Next take case (i). In this case the function f(z?) is of genus 47, order 
4i, and so is flat in a sequence of circles 
|je—l,| < daz. 
or z?—[2| < daz" |z+-C,| = 2dAzZ**+1+4 O(AZ"). 
Writing z* = w, (2 = y, this becomes 
w—y,| < 2dly,|-“4+4 O(\y,|-*. 
This region contains the circle 
w—y,| <dly,|-*+; 
and f(w) is therefore flat in this circle. 
The same argument disposes of cases (iii), (iv), (ii). 
6. Let f(z) be a function of finite order p, and let o< 1—43p. Write 
o = 1—}p’ so that p< p’, and take an integer j such that 
: 2 
air garg 
p’—p 


Then there is an integer 7 such that 


21-+2 
j 
Thus f(z’) is of order pj, and so, by Theorem B, is flat in a sequence 


Also, 


of circles z - 
z—€,| <da;*. 


Making the transformation z/ = w, we deduce that (2) is satisfied 
om 3 
in a sequence of circles 
2—t.| < d\u-1-ai 


and the index of A, is greater than o. This completes the proof of 


Theorem C, 


7. It is perhaps worth while to add that if y is any limit-point of 
the arguments of the points ¢,, the radius 6=y has the property 
that any sector containing it contains an infinite set of the points €.. 
By Bolzano’s theorem there must be at least one such radius. This 
is true both for Theorem B and for Theorem C. 








FURTHER STUDIES OF EMDEN’S AND SIMILAR % 
DIFFERENTIAL EQUATIONS 


By R. H. FOWLER (Cambridge). 
[Received 1 June 1931] 


1. Introduction. In a recent paper,* with perhaps some astro- 
physical applications, I have discussed in some detail the nature of 
the solutions of Emden’s equation 

agl?s +o" =0 (n > 1) (1) 
as € > 0, or of the equation 

0” +270" = 0 (n> 1) (2) 

as x — 00, which is equivalent to (1) when o = —4. This work arose 
out of my recent recognition, thanks to Professor Milne, of the astro- 
physical interest of continuing some much earlier incomplete work 
of my own on the nature of the solutions of the similar equationt 


- (w?0’) +-ba70" = 0 (3) 


as «> 0. Only integral values of n are contemplated in I, but most 
of the arguments transfer immediately to rational values. The num- 
ber of different cases presented by these equations in which the 
possible behaviour of the solutions differs and for which rather 
different arguments must be used is quite considerable, and the work 
quoted has not succeeded in giving an exhaustive analysis of the 
possibilities of every case. I attempt, I believe with success, in this 
third paper to fill in the main gaps, concluding it by repeating the 
complete results of all three papers in tabular form similar to that 
adopted in I. I do not seriously care whether these further investiga- 
tions can be defended as being of physical or astrophysical interest; 


a resolution of the grouping and behaviour of the solutions of these 
simple non-linear equations may perhaps be found of sufficient in- 
trinsic interest. But it is in fact likely that a detailed study such as 
this will prove of value in more than one physical problem as time 


goes on and the more difficult problems, which lead to non-linear 


* R.H. Fowler, M.N.R.A.S., vol. 91 (1930), p. 63. Quoted as IT. 
+ R.H. Fowler, Quart. J. of Math., vol. 45 (1914), p. 289, with references to 
| earlier work. Quoted as I. 


Ss? 
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equations of this type, come more into the foreground. Such a 
problem, for example, as the distribution of ions in and around an 
electric arc, allowing for recombination, seems to lead to an equation 
of this non-linear type.* 

2. Reduction of equation (3) to simpler forms. The earlier 
work in I refers to equation (3), but it is unnecessary to discuss an 
equation quite so general. (i) Suppose p > 1 and let 

dx/x? = dq/q’, q = (p—1)aP-, 0 = y/q. 
Then equation (3) becomes 

y"+bq7y"=0 (o' = (a+ p)/(p—1)—n—3), 
and q-—> 0 when x oo. (ii) Suppose p < 1, and let 


«., dq, q = (1—p)a?-?. 
aP 


Then equation (3) becomes 
6"+-bq7"'6"=0 = (0” = (o+p)/(1—p)), (5) 


and g > oo when x > oo. The constants b in (4) and (5) have the same 


sign as the b in (3). (iii) Finally, if p= 1, let g = logx. Then equation 
(3) becomes 6”+beMO"—0 (A=o+1). (6) 


In order to discuss all cases covered by equation (3) it is therefore 
sufficient to discuss equation (2) together with the equation 


6” + erg” — 0 (n > 1). (7) 


We shall here be concerned only with equation (2). The omission of 
the factor 6 in (2) and (7) can be arranged for by a simple change 
of scale of @, and is therefore trivial, but the sign of b matters. It is 
therefore necessary to include for completeness a discussion of 


0” —x70" = 0 
9” —eArgn - () 


* For studies of these equations reference must also be made to Emden’s 
Gaskugeln, especially Chap. XII where Emden discusses the general solution 
of his equation, not merely the special solution associated with his name. His 
results, while not, I think, complete, contain a fair proportion of the theorems 
given here, obtained in a rather different form. Where his results and the 
results of this paper overlap, they agree as they should. 

Mention should also be made of the work of Hopf, M.N.R.A.S. vol. 91 
(1931), p 653, in which he proves the theorems of IT and many of those of I by 
reducing the problem to the study of an equivalent equation of the first order. 
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We shall not discuss (9). Almost all these equations have already 
been adequately discussed in I when n = 1, and this case will not be 
touched on here. 

In all cases, what is required is a full description of the possible 
behaviour of real continuous solutions of these equations, possessing 
continuous first differential coefficients as x > oo. Such solutions will 
be called proper solutions. It will be assumed that n = p/q, where 
p, 7 are integers relatively prime. The parity of p and q is important 
in so far as it controls the behaviour of 6” as @ crosses or tries to 
cross 6 = 0. If p, g are both odd we shall be able to take (—@)" = —@” 
and to say that n is odd. If p is even and q is odd, (—6@)" = 6", and 
we say that n is even. If p is odd and q is even, (—@)”" cannot be real 
and proper solutions cannot ever cross # = 0, still less oscillate about 
60. This case is of rather less interest (mathematically). The chief 
interest lies in the equations (2) and (7) rather than (8) and (9) and 
in n odd rather than in m even, for it is only for (2) and (7) and n odd 
that oscillating solutions can occur. 

3. Outline of known results. Detailed statements of the be- 
haviour of the solutions in all cases will be given at the end of this 
paper. It will perhaps be best, however, to start with a simplified 
description of known results for equation (2). It has been shown in 
[ and II that the nature of the solutions of equation (2) is largely 
controlled by the signs of the two expressions o+-2 and o+n-+1. If 
o-+-2=0 or o+n+1 = 0, equation (2) can be reduced to | 

d*@ «dé 
dt? dt 
d*y dy 
di?" dt 


+6" =0 (t= log), (10) 


+y"=0 (!=log2, y="0), (11) 


respectively. The first equation has been fully discussed in I, see 
Table VI there; all its solutions must oscillate as ¢t > 00, with oscilla- 
tions of a specified form. The second has been fully discussed in IT; 
all its solutions are ultimately monotonic as ¢ > 00, again with speci- 
fied forms, and specified relative arrangements. There remain the 
three cases of inequality: 

(i) If o+2> 0, then o+n+1> 0a fortiori. This case is fully dis- 
cussed in I. We may summarize by saying that the behaviour is 
similar to that of o+2= 0. 

(ii) If o+-n+1<0, then o+2<0 a fortiori. This is a case of 
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astrophysical interest, fully discussed in II. The behaviour is similar 
to that of o+-n+1=0. 

(iii) c+ 2<0,0+n+1>0. This is the most intricate and mathe- 


matically most interesting case. It is our main subject of study here. 


It subdivides into three according as 20+n+3<,=,>0. It may 
sometimes possess simultaneously solutions which oscillate and others 
which are monotonic and tend to zero as x > 00. 

4. The equation 

6"+-2°9"=0 (o+n+1>0,04+2<0,n>1,nodd); (A) 
various transformations and notation. Certain equivalent forms 
of equation (A) are here of greater assistance. On substituting 

= Ax™ (A, @ constant), 
we find that there are two special solutions of the form 
, - 
CO+ 4 
(~ — ). (12) 
n—1l, 
These are proper solutions since 0 < @a< 1. We therefore write 
1 
=[aw(1—a)|"" ‘xv, (13) 
and form the differential equation for v, which is 
xv" + 220’ —aw(l—w)(v—v") = 0. (14) 
We can further eliminate x from (14) by writing t=logx, when 
(14) becomes 
vp” + (2a—1)v’—ev(1—v"-") = 0 (c= w—w’). (15) 
The special solutions are now v= +1, and, of course, v=0. No 
other solution can touch these, for it would then be identical with 
the solution touched. 

In (15) dashes denote, of course, differentiation with respect to f, 
and we shall always use them with this meaning relative to which- 
ever independent variable is in use at the moment. We shall retain 
throughout the standardized relationships 

= = 1/f (t,x > oo, €> 0), (16) 
. 1 1 

6 = ay = ey =[w(1—w) |" 'a7v = [o(1—w) |" *e™. (17) 

The sign of 2a~—1 changes when mw = }, i.e. when 

2o-+n+3 = 0. (18) 


The astrophysical example is o = —4, n= 5, first investigated by 
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Schuster. The discussion of this case is elementary and the results 
more or less known, but a short survey from the present viewpoint 
will serve to introduce the general cases 

2a—1> 0, 20+n+3 < 0, 
20—1< 0, 20+n+3> 0. 
5. Equation (A). The special case 20+n+3=0. We have 
to discuss the (v,t) equation 
v” = }(v—v") (t > 00). (19) 
This integrates to give 
1/2 — 


Dy 


Ly? v1 D. (20) 


ca _ yh 
4(n+-1) 
Since n > 1, v-> + 00 implies v’? + —oo, and so is impossible. We 


can write ie 





= dt ° 
Eppes l 1) 
and v can at most oscillate to and fro between the greatest and least 
roots of 1 

aad ila 22) 

When D = 0, Pe 





= —hdt, 


» l — vt) 
n+l 


where the ambiguous sign has been chosen so that t-> oo. With the 


help of the substitution 
2 
v"-1 = sin*A, 
n+1 

this is easily integrated in the form 

9 2p—4(n-1t 

Bg an ee 3? (24) 

n+1 {1+ Cre-tn D2 


C being a constant of integration. The corresponding curves in other 


coordinates are = 

wD) Y2y4(n—1 - 

2(n+ 1)Craie |" (25) 
fain _4 O2}2 ? 


- fl + C2 eine 


Hn 1)? Jn 
{1 + aT f 


(n-1)C2 i 
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These solutions are finite at the origin and are therefore of Emden’s 
type. They were, of course, first given by Schuster for n = 5. They 
are monotonic everywhere, decreasing as € increases, and tend to 
zero like A/E& as E> «wo, x > 0. 

The (v, ) Schuster-Emden curves form a one-parameter family with 


an asymptote easily shown to be 
1 


+v=(* sl —) 


The exceptional solutions v = + 1 therefore lie in the region covered 
by the Schuster-Emden solutions, which in (v,/) coordinates are all 
the same curve with different origins for ¢. Since, as is easily proved, 
y(t) > 0 for t <t and v’(t)<0 for t>t 
Emden curves intersect once and once only. They therefore form 


men max? @ pair of Schuster- 


a grid with the necessary properties for use in analysing the other 
solutions like the Emden-curves of the cases treated in I. 

The properties of these other solutions can of course be directly 
derived in detail from the equation (21), but it is instructive to see 
how far we can proceed with the grid. For this purpose one must 
know something about what the other solutions do as t > o, and it 


is obvious from (21) that they oscillate about at least one of v= +1, 


0, —1 and cannot tend to zero, infinity, or a finite limit. We argue 
here in (0, €) coordinates in which the diagrams are easiest to com- 
prehend. 

On referring to Fig. 1, let P, be any point below the envelope 

1 

6= {h(n+1)}""é-?. Then by using, as in II, the argument that a 
solution which is not an Emden curve may never touch any one of 
the family of Emden curves, we can prove that any solution such as 
P,1 which starts inwards above the upper Emden curve through P, 
must pass above the envelope. Any solution such as P,2 which starts 
inwards between the two Emden curves through P, must lie as shown 
in Fig. 1 and tend steadily to infinity as € > 0, always lying below 
the envelope. Lastly, any solution which starts inwards below the 
lower Emden curve through P, must cut the €-axis at-some point 
between 0 and the abscissa of P,. It cannot cut the Emden curve 
itself without touching an Emden curve and cannot lie continually 
between € = 0 and the Emden curve without making v(t) > 0. 

Again, any solution which passes above or starts above the en- 
velope may not stay above the envelope, for v(t) >k>1 as t+ o, 
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and therefore must cut the envelope again, and so the Emden curves, 
in such a way that it must cut the é-axis again as in the last case. 


Envelope 














\ 
Fic. 1 


It can never touch the envelope. Further, solutions which cut the 
€-axis at no matter what angle, whether from above or below, must 
lie near the axis relative to the Emden curves there as shown at P, 
and must pass above the envelope again. For obviously as P tends 
to the €-axis the slope of the two Emden curves through P both tend 
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to zero, the Emden curves of this equation having nowhere a zero. 


We have thus proved: 


THEOREM I. As & - 0, solutions of the equation (A) (2o-+-n+-3 = 0) 
other than Emden’s solutions must either (1) oscillate infinitely and 
infinitely often, cutting infinitely often both the curves 

1 
g= (le, (28) 
8 / 

or (2) tend monotonically to +00 (or — oo) always lying below (or above) 
the relevant curve (28). The slopes of the Emden curves through any 
point P between the curves (28) enclose between them the slopes of all 
the solutions which tend monotonically to +00. Solutions through all 
other points or of any other slope oscillate infinitely as € + 0. 

To prove more than this one must of course go back to equation 
(21). It is then simple to prove by studying the roots of the expres- 
sion under the square root that the oscillating 9-solutions are solu- 
tions in (v,t) which oscillate continually and with constant form 
about all three of the lines v= 0, +1. These solutions therefore 


satisfy the asymptotic form 


Cag MEF,» fA Winds. '), (29) 


3) 


max 


é,,, = Bé, (B <1), (30) 


where &,,,, €, are any two consecutive zeros. B is a definite function 


of A but of a form not explicitly determinable. The monotonic 
6-solutions are solutions in (v,t) which oscillate continually and with 
constant form about the line v = +1 (or v= —1) only. These mono- 
tonic (!) solutions therefore satisfy the asymptotic conditions 

1 


+ ]\n 
CyEA<OE)< OE (C,>0,0,<(2)""), 


C, and C, are functionally related, and 6(€) takes repeatedly all 
1 


(31) 


P ; n-+-1\n-1 
values in the range as €>0. C,—>0 when C,—> ) , and 
. “ 8 
1 = 
C, > (4)"" when C,- (})"-1. The limiting forms of these solutions 


are the special solutions v(t) = +1 i.e. 


1 
= +(e. 
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6. Equation (A) (20+-+3< 0). We can discuss this equation in 


the form vy” +bv’+ce(v" —v) = 0, (33) 


where 6=2H~—1>0, c=mw(l—w)>0, and too. Except for 
certain details arising from the extra complication of the function 
v"—v in place of v” the discussion follows very closely the discussion 
of y"+y’+y" = 0 in II, and proofs which are mere transcriptions 
of proofs there given will not be repeated. 

Consider first the limited case in which it is assumed that there 
exists no fy such that v(t) << 1 fort >t). For this case we can prove 


THEOREM II. If there exists no value ty such that \v(t)| <1 when 
i >to, any proper solution v = v(t) of equation (A) (20+n-+3 < 0) ts 
such that as t-> o either v(t)» +1 or v(t) > —1 with or without an 
infinity of intersections with v = +1 or v= —1 respectively. 

There are here four a priori possibilities of behaviour relative to 
the line vy = +1, with corresponding ones for v = —1: 

(i) There exists a t, such that, when ¢ >t), v(t) > 1+k, k>0. 

(ii) For any positive k there exists a ¢, such that, when ¢ > fy, 
v(t) > 1 and v= v(t) cuts v= 1+k in a {t,}.* 

(iii) v = v(t) oscillates finitely or infinitely about v = 1 cutting both 

1+k and v = 1 in {t,}’s. 

(iv) v(t) > 1 as t-> oo with or without cutting v = 1 in a {t,}. 

The theorem consists in the impossibility of sub-cases (i), (ii), (ili). 
Sub-cases (i), (ii) are both impossible by proofs as in II which will 
be omitted. Sub-case (iii) is also impossible, the proof in IT requiring 
some amplifications. Assuming for the moment that between ¢, and 
t,.;, any pair of consecutive intersections of v = v(t) and v = 1, there 
is just one value of t, 7,, at which v’(r,) = 0, then it follows just as 


in IT that [v’(t,) |? > f{v(z,)} > [v'(t,.4))?. 


1 
with the same essential properties makes no difference. It follows 
further, just as in II, that |v’(¢,)| +1 + 0 is impossible, and therefore 
se [oO (34) 
and that (34) implies that v%(r,) > 1. (35) 


2 
The substitution of a more complicated function for wey a 
n 


* By {t,} we denote an enumerable sequence of values of ¢ tending steadily 
to infinity. 
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Let us return next to the uniqueness of 7,. This presents no diffi- 
culty on ares on which v(t) > 1. On ares for which v(t) < 1 either of 


the configurations shown in Figs. 2 and 3 are possible. It is further 
easy to show, by considering the sign of v’(t) near a zero of v’(t), that 


Nis Ts +1 — Ats+2 





v=l 























Fic. 3 
the only possible way in which v(t) can vanish more than once 
between ¢,,, and ¢,,, is as shown in Fig. 4. The signs of v”(¢) in the 


—fXse0 


different regions are indicated in the figures. 





cae \ts+i a 
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But then we can conclude by the same argument, which only 
requires that v’(¢) should preserve its sign on the arc used, that 


[etd P >for}, — Avlrs)} > [v'(tss2) P- 
But now, on applying the same argument to the intersections of 
» = v(t) with v = —1, we can conclude that 
S(t) > [e'Had P > foliad} > [ed P > Kors sd} 
the same function f(v) being involved. The essential conclusions, that 
Iv'(tesa)| > le'(tesa)|, 0°.) > 0, 
remain valid. 

Now that we know that v’(t,) > 0, we can return and re-examine 
whether ares such as those shown in Figs. 3, 4 can continue to survive 
as t-> 00. Suppose they do survive. Consider the arc starting from 
t.., in Fig. 3, on which, for large enough s, wv’(t,,,) like v’(t,) can be 
made as small as we please. On the arc ¢,,,t (¢ < t®,,) 

—v" = br'+c(v"—?). 
Both terms on the right are positive, since v(t) > 0 and —1 < v(t) <0. 


Therefore 
—_ vy” he —c(v—v"). 


All along this are v = v(t) turns over more quickly at corresponding 


heights than the solution v = V(t) of 
—v" = —c(v—v"), 

which has the same starting-point and starting-slope. The curve 
v = v(t) must therefore attain its maximum at which wv’(¢) = 0 nearer* 
to v= —1 than v= V(t). But by direct integration 
i 7 n+1 7 2 aa l 
Vea) 20(LPOP" WP) ot | _ 

| n+l 2 '2n+1)J 
and the root V(r) of this equation nearest to V = —1 approaches 
V = —1 as a limit, as V'(t,,,) + 0. The first vanishing-point of v’(#) 
must lie even nearer. Hence the forms shown in Figs. 3, 4 are im- 
possible for t >t), and v = v(t) reduces to a curve such that either 
v(7,) > +1 or v(7,) > —1. Since, by definition, v(7,) represents the 
points of maximum excursion from v = +1 or v = —1, we can con- 
clude that on these solutions either v(t) > +1 or v(t) —l ast o. 


LV (7)}, 


* It does not appear to be simple to prove this statement rigorously except 
by a reductio ad absurdum, starting from the hypothesis that the zero of v’(¢) 
lies further from v = —1 than the zero of V(t). 
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A solution may start (as we shall later prove rigorously) by oscillating 
about all of » = 0, +1, but, if it continues ever to lie above v= +1 
or below v -las t> oo, then v(t) +1 or v(t) > —1 respectively. 
This completes the proof of Theorem IT. 

One further point should be made, and that is that configurations 
such as those shown in Fig. 4 are everywhere impossible. This is most 
simply proved by showing that, to the left of an arc such as ?8,,7,,, 
the next arc in the region —1 < v < 0 is everywhere steeper at corre- 
sponding heights and therefore must again intersect v = 0. The proof 
follows at once from the usual inequalities based on the sign of b and 
can be left to the reader. 

The remaining alternative is that there exists a ¢, such that 

v(t)|} <1 for t >t). To cover this we prove 

THeoreM III. If there exists a t, such that \v(t)) <1, when t > to, 
then any proper solution v = v(t) of equation (A) (20-+-n+3 < 0) ts such 
that, as t—> «, either v(t) > +1—0, v(t) > 0, or v(t) > —1+0 mono- 
tonically. 

We can observe at once that proper solutions restricted to |v(t)| < 1 
cannot cut v= k (0<k <1) twice and lie between v =k and v= 1 
between the intersections; nor cut v —k (0< k <1) twice and lie 
between v -k and v= —1 between the intersections. It follows 
at once that for any given k (—1 <k < 1) there exists a ¢, such that, 
when ¢ >t), either v(t) > k or v(t) <k, and therefore that there must 
exist an / (— 1 <1 < 1) such that v(t) > / monotonically. Proper solu- 
tions for which /—0 will be found to exist (Emden’s solutions). 
Proper solutions for which /= +1 may also exist, but they are, in 
no essential, different from those allowed by Theorem I. Proper 
solutions for other values of / are impossible. For then 


vw” +bv’ ~ a(~ 0), vy’ +bv ~ at, 
x x 
pelt ~ tel”. v~ t, 
b b 


which contradicts v(t) +1. This completes the proof of Theorem III. 

When v(t) > 0 monotonically the form of the solutions, if they 

exist, can be determined at once by Hardy’s theorem.* For, on 
writing v’ = 7, the equation becomes 

dy 

Che an 

dv 


Hardy, Proc. London Math. Soc., vol. 10 (1911), p. 451. See also IT, p. 69. 


+ by—c(v—v") = 0, (36) 
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and the form of solutions is required for which v > 0, » > 0, t > o. 
It is easily found that such solutions must be of the form 

v(t) ~ Ae-™, 
where A is a constant of integration. Such solutions satisfy also 
OE) > C. 


They are Emden’s solutions, and their actual existence and unique- 
ness may be established as in II. We embody these results in 

THEOREM IV. Emden’s solutions of equation (A) (20+n+-3 < 0) 
satisfying , 

lid Ag)>+C (E+) 
where C is arbitrary, exist and are unique. If v(t) > 0 as t > co on any 
proper solution, then must 
v(t) ~ Ae-* 
and the solution must be an Emden solution. 

A more exact study of the solutions v(t) > +1 is also possible and 
will be given later. It is best to discuss first equation (A) when 
2a—1 <0, t+ , for, since it is then the same equation as when 
2a—1 > 0, t > —oo, knowledge of one case can be used to help with 
the other. 

7. Equation (A) (20+-n+3 > 0). The form of the equation is now 

v” —bv’ +c(v"—v) = 0, (37) 
where b = —(2a—1) > 0, c= a(1—w) > 0, and t+. We prove 
first 

THEOREM V. Equation (A) (2¢o-+-n+3 > 0) possesses no proper solu- 
tions, other than v(t) = +1, such that either v(t) >1 when t >to, or 
v(t) < —1 when t > ty, or v(t) + 1—0, or v(t) +> —1+-0 monotonically. 

Solutions obeying the restriction v(t) > 1 must be ultimately mono- 
tonic and their asymptotic behaviour can therefore be studied by 
Hardy’s theorem applied to 


ng —by+c(v"—v) = 0. (38) 
av 


It is easily shown by examination of each alternative that v(t) > 00, 
v(t) > 1, and v(t) > 1+0 all lead to contradictions. The details are left 
to the reader. We can also conclude by the same argument that 
v(t) > 1—0 monotonically is equally impossible and also the alter- 
natives concerning v = —1. Hence the theorem. 
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We prove next 

THEOREM VI. The only proper solutions, other than v(t) = +1, of 
equation (A) (20+n+3 > 0), 
which satisfy the restriction that 

—l<v(t)<l 
when t >t), are Emden’s solutions which exist, are unique, and satisfy 
v(t) ~ Ae~™, AE) > C (€ > 0). 

We know already that v(t) > +1, —1 are impossible. Two con- 

secutive zeros of v’(t), if they exist, must lie on opposite sides of 


| 











Fic. 5 


v = 0, and therefore must be arranged as in Fig. 5. But v” > 0 at P 
and v” < 0 at Q, so that this is impossible. Hence all such v = v(t) 
are ultimately monotonic and therefore v(t)>/1 (—1<l< 1). If 
v(t) +1( 0), then 

v” —bv’ ~ a( 0), v’ —bv ~ al, 
xt 


—bt 
-€ . 
b 


d 


54 (0 bt) — ote, ve%t_C,~C, 
C 


which contradicts v(t) > 1. 
The only alternative is »(t) > 0 monotonically. Just as before, we 
can prove that Emden’s solutions exist and are unique, and that, 


when v(t) > 0, v = v(t) must be one of Emden’s solutions. This com- 


pletes the proof of Theorem VI. 

The remaining possibility is that v= v(t) cuts v=1 and (or) 
v = —1 infinitely often in a {t,}. In this case we shall prove that 
- x(t) must ultimately oscillate infinitely and infinitely often, cutting 


9 = 
0, +1 in {t,}’s. Combining this result with the results 


all three of v 
of Theorems V and VI, we thus prove 
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TueEorEM VII. If a proper solution of equation (A) is not one of 
Emden’s solutions or one of the special solutions v(t) = +1, then, as 
t > 0, v = v(t) oscillates infinitely and infinitely often, cutting v = 0, +1 
in {t,}’s. 

The arguments to be used are exactly the same as those of the 
corresponding part of § 6, except that the sign of b is changed so that 
the sense of every inequality is reversed. We can therefore prove at 
once that, if {t,} is the sequence of intersections of v = v(t) with v = 1, 
walk 1)< Io's) 
and further that, as t, > 00, 

|v’ (t,)| > 00. 

We have already proved in § 6 that oscillations about v = +1 (or 
v= —1) and v= 0, +1 cannot alternate; the same argument shows 
that any oscillation which succeeds an oscillation about v= 0, +1 
as too is also an oscillation about v= 0, +1. It therefore only 
remains to show that since |v’(¢,)| > 00 it is impossible for the oscilla- 
tions to remain confined to v=-+1 (or v=—1). Consider an 
oscillation below v= +1 starting with a given sufficiently large 
value of |v’(t,)|. On this are 

v” = br’+c(v—v") 
and bv’ <0, so that v" <ec(v—v"); 


the negative slope diminishes less rapidly than on the curve v = V(t) 
which is the solution of 


v” = c(v—v") 


with the same starting-point and starting-slope on v= 1. Therefore 
at corresponding heights on these arcs 
\v'(t)| > |V'"@)|. 

But by direct integration we can show that v= V(t) must inter- 
sect v=0 when |v’(t,)| is sufficiently large. Therefore v= v(t) 
must also intersect v= 0, and therefore v= 0, +1, for sufficiently 
large |v'(t,)| and therefore when ¢t >¢t). This completes the proof of 
Theorem VII. 

8. Properties of the Emden solutions. Consider first the case 
2c+n+3 <0. Then the Emden solutions in (v,¢) coordinates are 
those solutions which satisfy v(t) ~ Ae-™ as too. They are ob- 


viously all the same curve with variable origins for ¢. We require to 
3695.2 
T 
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know how they behave as we follow them backwards towards 
t- —oo. In this direction they are controlled by Theorems V—VII 
and therefore, as t-> —oo, v= v,(f) either satisfies v,(t)> 0 and is 
again an Emden solution, or else oscillates continually about 
v=0, +1. In order to turn over so that v,(t) can start diminishing 
as we trace v = v,(t) backwards it is necessary for v = v,(t) to rise 
above v= 1. [It is tacitly assumed that v,(¢) is positive for large t. | 
Now it is impossible for v = v,(t) to cut v= 1 more than twice, so 
making a loop below it before it cuts v= 0. For on the are which 
starts off below v = 1 the slope must be everywhere steeper at equal 
heights than on the Emden branch ¢ > +00, and this is impossible 
if v;,(t) vanishes for positive v. Thus the only remaining possibilities 
are those shown in Fig. 6. And in fact it is easily seen that the dotted 
alternative is impossible. For in that case, on the arc ty, +00, 


v” = —bv’+c(v—v"), 


v” >c(v—v"), 
—v'v" > —c(v—v")r’" 
1 
$v’ (ty) |? > [ (v—v") dv. 
0 


Similarly on the are —oot, 


But from the arc t,t, |v’ (to) | < |v’ (é,)|, 

and these are contradictory. Hence the case of the dotted curve is 
impossible and the full curve of Fig. 6 shows the nature of the Emden 
solutions. We observe that, between the last zero z and o, v’(t) 
vanishes just once, and therefore that any two curves of the family 
of Emden curves stopped off at their first zero intersect once only, 
and that the family has a single envelope v=a>1. In this case 
therefore the family of Emden curves still forms a grid which can be 
used to define the behaviour of the other solutions. 

When 2c+2+3>0 similar arguments may be applied. The 
Emden curve v = v,(t), when continued back towards t > —oo, must 
either be again an Emden curve or tend to v= +1 or v= —1 with 
or without oscillations. The argument just given shows that it cannot 
be an Emden curve again. [This happens only when 20+-n+3 = 0.] 
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It cannot approach v = —1 without cutting v = 0, and such an arc 
is impossible because its slope must be less than the slope of the 


AN Pe Be vzl 


“, 
Af to 
\ TS as vip 
—_ = v=0 
\ ei 
v=-l 
a 5 ae so 


——> 00 E->0 
Fic. 6. Emden solutions for 2o6+n+3 < 0. 














Emden are at corresponding heights. Its nature is therefore fully 
shown in Fig. 7. 

Emden’s solutions never vanish. It is at once evident that such 
curves may have an infinity of envelopes and any pair an infinity of 
intersections before they reach v = 0 = 0, so that they no longer have 
the properties necessary for use as a grid. And in fact the other 
solutions, except the special ones v= +1, are now all of one type 
as > +00. 


eM | 
0 











t —>-© t—>+ 00 
iis: ——>0 


Fic. 7. Emden solutions for 20+n+3>0. 


We can now use the grid of Emden solutions, but only when 
2o+n+3 <0, just as we have already used it for 20+n+3 = 0 and 
for co+-n+1<0, o+ 2<0 in II. As a result we can establish the 
final points required for the following theorems, which also recapitu- 
late information given previously. 

THEorEM VIII. Jf n= p/q, where p, q are odd and if 

o+n+1>0, 20+n+3< 0, o+2<0, 
then no proper solution of equation (A) 
0” +279" = 0 (a —> 00) 
T2 
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or of the ny 


agg) ter te =9 E50 
can continue to oscillate about 6 = 0 
Such solutions must take one or other of the asymptotic forms 
6= C— — (a ~ 1), [ E-solutions] 
(o+1)(o+2) 
where C is a constant of integration, or 


Lo 


1 
a~ +[— eee | | M-solutions| 
n—1)2 


of which the former are exceptional. There are two special M-solutions, 


1 Lé 
9 — tb. —(e+2)(o+ n-+-1) ie 8 
7 (n—1)? ; 


THEOREMIX. Through any point P, 0,, €, (0, > 0) below the envelope 
of the E-solutions of equation (A) (20-+-n+-3 < 0) namely 


o+2 1 ; 
Dns mg ™ (m = j= +2 ae 


(n—1)* 


pass two E-solutions with slopes* a, B, a>. (i) If «>6(&) >, 
i= al is an M-solution, monotonic when & < &,, and 0(€) > +2 as 
E€>0. (ii as 0’'(€,) =a, or 0'(€,) =f, the solution is an E-solution. 
(iii) If ba (€,) >a, or 0’(E,) < B, the solution cuts 0 = 0 again at a point 
&, such that 0 < & < &,. 

Any proper solution through a point P, 0,, €, above or on the envelope 
cuts 0 = 0 again at €, where 0 < &, << &,, except that through a point P 
on the envelope there passes the single E-solution which touches the 
envelope there. 

Similar relations hold for negative 6,. 


9. More exact asymptotic forms when 2c+n+3<0. Let us 
consider more closely the solutions for which v(t)» +1 in this case. 
Writing v = 1+v the equation becomes 


v" +b’ + (n—1)ev[1+f(v)] = 0, (39) 


* In the enunciation of this theorem slopes are reckoned positive when 
6(é) increases as € > 0. 
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where f(v) is a rational function of v and v” such that f(v) > 0 as v > 0. 
It is almost obvious, since we know already that v(t) 0 on the 
solutions in which we are interested, that the asymptotic forms of 
these solutions must be the solutions of the equation 


v” +b’ +(n—1)ev= 0. (40) 


A rigorous formal proof, however, seems to require the methods used 
in I for discussing standard form II. The arguments there used can 
be extended at once to cover what is required here, and we may 
therefore proceed to use the solutions of (40) as the required asymp- 
totic forms.* 

Solutions of (40) are of the form 


_ —b4%b2—4(n—1)c 
7” 2 PS ig 2 —® 
Thus v = v(f) is ultimately monotonic if b? > 4(n—1)e or 


(2o+-n-+43)? > —4(n—1)(o+2)(0+-n+1). (42) 


where r (41) 





In this case 


v(t)—1 se v(t) — eh(20+n +3M elu +e + eet eat] (e,, €—> 0), (43) 


where [(20-++-n-+3)?+4(n—1)(o+2)(o+-n-+1)]! (44) 


= — 9 


2(n—1) 





and the second term in (43) may only be retained if the first vanishes 
identically. If, however, 

(20-+-n+3)? < —4(n—1)(0+2)(0-+n+1), (45) 
then v =v(t) oscillates continually about v= 0. If ¢,, t,,, are con- 
secutive zeros and 7, the abscissa of the maximum, then 

v(7,)—1 =v(7,) ~ Cet@otn+sy. (46) 
where C is a constant of integration, and 
iene 4n(n—1) 
| [—4(n—1)(o + 2)(o-+-n+1)—(20-+n+3)] 





* An alternative line of argument is to construct by an iterative process 
exact solutions of equation (39) of the required form by using the solutions 
of (40) as starting-points and showing that these solutions of (40) are then 
necessarily the asymptotic forms required for the solutions of (39). This is 
straightforward and contains no point of novelty. 
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10. More exact asymptotic forms when 20+n7+3>0. The 
necessary theorems to detail to us the forms of the infinite oscillations 
of the solutions other than Emden’s and the special solutions in this 
case have already been given implicitly in I. Theorem VI (a) of I can 
be repeated for equation (A) in the following form: 
If {x,} be the sequence of zeros and X,, 0(X,) the sequence of maxima 
and minima on any oscillating solution, and if x,,,/%,—> 1, then 
JAX =X, "(e+ 0), (48) 
418 = ty 1d (« > 0). (49) 

It is not, however, proved in I that if a+ 2 < 0, 20+n-+3 > 0, then 
s+1 >a ? 

To supply the missing step it seems to be simplest to extend the 
results for the standard form II of I, Theorem VI (b). This theorem 
states inter alia that 

If {t,} be the sequence of zeros and 7,, v(7,) the sequence of maxima 
and minima on any oscillating solution of 


v" —bv’ cv" = 0, (50) 


x / x 


s 


and if t,,,—t, > 0, then 


u(r.) = exp| (=5-+«)¢] (e > 0), 


\\n 


teat, = exp| (—S5>+<) | (> 0). 


\\ +3 
In particular, if as here b > 0, then t,,,—t, > 0 and relations (51), (52) 
hold good. 

It will be observed that (50) is not the desired equation and we 
must repeat the arguments with v”—v in place of v", but we know, 
to start with, that v= v(t) oscillates about all three of v= 0, +1, 
having maxima only above v==1 and minima below v= —1, and 
that |v(7,)| > co. We can confine attention to b>0,¢c>0. 

The correct v-equation may be written in the form 


deat atest —oy. 
0 


If t) be any point on the arc ¢,7,, then for fj <t <r, 


“s ~2U ~2bt 
e- 207s < e-2bt < e-2lo, 
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and therefore 


2ce-27e(y" —v)v’ < — . (e-%y’)? < 2ce—2lo(y — yy’. 


vs(Ts) 
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Integrating this over the range tyr, we obtain 


Qce-WTs {[v(zs)]}"* a [v(r,)}? _ [e(to)]"* re eto 
a 








n+1 2 n+1 2 
< e-2%hif y(t.) 2? << 2ee-2% isco wit [ora r et Pe evo 


These inequalities hold for all values of t, in the interval and therefore 


2Qce-2(rs—) [o(7,)]"** = [v(7,)]? vs gntl =a 











n+1 2 m+1' 2 
, ' [v(7,)]” +1 [v(z,) |? prt 2 
<(v')? <2e a” 8 +5}: (54) 


Let us now compare the v = v(t) curve with the curve 


v=folp), (55) 
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where f,(#) is the perfectly definite real function defined by 
n+1 2 n-+ 
(or _@ £0" | fe\ty, ay 


\n+l 2 n+l 2] 


taken with a suitable origin. Let 


| 
J\n+1 2 J 
0 
then fg{K(Q)} = Q. As t goes from ¢, to 7,, % goes from 0 to K(Q) 
and v(t) goes from 0 to v(7,), with a similar relationship for the second 
part of the wave. On comparing the inequalities (54) with (55) and 
(56), we obtain 


+18) fg = K(Q): (57) 


i Qn +1 Q? Be tl 


Qce-2rs N— (pb )P< 2c; 


we take the square root of these inequalities (all square roots being 
positive) and integrate from ¢, to 7,, thus obtaining 


PP 1 —e-We-0)} < K(Q)< (20)!(r,—ty) (58) 
Similar arguments applied to the arc 7,,,, show that 
j r (2¢)* na 5 
(2c)¥(t,41—1,) << K(Q) < {eer — J}. (59) 
) 
It is immediately obvious that we have only to prove that K(Q) + 0 


as Q(= |v(r,)|) > 00, in order to deduce that ¢,,,—?, > 0 and so com- 
plete the investigation of this section. By (57) 


1 
|. 1) ae —4 
K(Q) oe (n+-1) | f{y—ans_”™ I 1 1 (1—a)| dx, 


Qia-)) , \ 2 Q | 


1 
_(n+1)! emt, eA 1—2? )-! 
x Qir-D [a-s 5 i ras 2 Q” -1 at da 
0 


Over the range 0<x<1 the coefficient of 1/Q”-1 in the second 
factor in the integrand is bounded, and therefore, as Q > 0, 


K(Q) ~ ~ oe ” { (1—a+1)-* da + 0. 


0 
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These results may be embodied in the following: 
THEOREM X. If as x > © {x,} be the sequence of zeros and X,, 0(X,) 
the sequence of maxima and minima on any oscillating solution of 
equation (A) (20-+n+3 > 0), then 


\(X,)| =X, **° (e+), (48) 


20 


€ 
Sua, =, OF (e+ 0). (49) 


These expressions may also be obtained in the more exact forms* 


Co 


\(X,)| ~~ AX, **, (60) 


1 1 
r(5)r ( aa - 
1 y 2 n + 1 ~nt+3 61 
8 ~ sina] (seq) a . (61) 
r{-+——~ 
(5 n+ i) 
11. Tabulations of the behaviour of the proper solutions of 
the equation 
a"+27"=0 (n>1,n=p/q, p,q odd). (2) 
We have now a fairly complete knowledge of the behaviour of any 
proper solution of equation (2) which may be exhibited as follows: 
Table I 
The form as x - © of proper solutions of 
0” +276" = 0 (n> 1, n= p/q, p, g odd). 
Case I.¢ o+2>0, 20+n+3>0, o+n+1>0. 
All solutions must oscillate about 06 = 0 with the asymptotic forms 


rch) f © 


~n+3 
r 1 1 x, : 
(3 sy, J 


where A is a constant of integration, {X,} the sequence of stationary 
points of 0(x), and {x,} its sequence of zeros. 


o 


\0(X,)| ~ AX, **°, 





1 2 \t 
X41, ~ Aiea ati 





* By the methods outlined in I, §§ 12-20. 
+ Sources: ¢ +2> 0, I, Table III; o+2 = 0, I, Table V, with a mistake cor- 
rected, which occurred in transcribing the correct results of I, Theorem VI (6). 
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Case II. o+ 2<0, 20+n+3>0, o+-n+1>0. 
There are solutions of three distinct types, some oscillating, some 
ultimately monotonic. 
(i) Special solutions (with no arbitrary constants) 
1 
O= + ea ”| ae. 
(n—1)? 
(ii) Emden’s solutions (with one arbitrary constant C) 
aCnget? em 0). 
(¢+-1)(o+2) 
Emden’s solutions never vanish for any value of x, and are everywhere 
monotonic. 
(iii) All other solutions (with two arbitrary constants) must oscillate, 
having oscillations of form (P), Case I. 
Case III. o+2 <0, 260+n+3=0, o+n+1>0. 
There are solutions of four distinct types, some oscillating, some 
ultimately monotonic. 
(i) Special solutions of form (Q), Case IT. 


(ii) Emden’s solutions (with one arbitrary constant C) 
1 


Fi (n+ 1)0° . ag 
{l + C2/ain ~~? 
The exact form (S) conforms to (R), Case I. 
(iii) Monotonic solutions (with two arbitrary constants) having no zero 
which tend to +-c00(—00). These solutions (+00) satisfy the relations 


C12! < W(x) < C,2%, 


@= C— 


1 T) 
- ? 1\n-1 ( 
1>0, O<("E)"",  G=NO%) 


/ 


and 6(2) takes all values repeatedly in this range as x —> oo. 
1 1 = 


1 
n+1)\n-1 -1 . n-1 
C,> 0, as Cy, > =) ;O,> (4)" ,as C,>(}) . 


(iv) Solutions (with two arbitrary constants) which continue to oscillate 
about zero as x -> 0, of the form 


\(X,)| = AXt, A >( vu) 


8 
%5i1 = Ba,, B>1, B=9(A). 
All these relations in Case III are exact, not asymptotic. 


1 
a 
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Emden’s solutions form a grid separating the monotonic from the 
oscillating solutions. 
Case IV. o+-2<0, 20+n+3< 0, o+-n+1>0. 
All the solutions are ultimately monotonic as x-> 00. They are of 
three types. 
(i) Special solutions of form (Q), Case II. 
(ii) Emden’s solutions (with one arbitrary constant) of form (R), 
Case II, but here these solutions have an infinity of zeros as x > 0. 
(iii) Ultimately monotonic solutions (two arbitrary constants) which 


satisfy 





1 
—(o+2)(o+n+1)]aa -2? 

a2) ~ +. |=. v) 

These solutions have an infinity of zeros as x > 0. 

Emden’s solutions form a grid separating off the solutions which have 
had their last zero as x > «© and become monotonic, from those which 
have more zeros to come. 

Case V. o+ 2< 0, 20++3<0, o+n+1=0. 

All the solutions are ultimately monotonic as x-> 0. They are of 
only two types. 

(i) Emden’s solutions (with one arbitrary constant) of form (R), 
having zeros as x > 0, as in Case IV. 

(ii) All other solutions (with two arbitrary constants) have the asymptotic 
form as x > 


1 ne 
ise ++ Tetie| iin 


where D is one of the constants of integration. These solutions all have 
an infinity of zeros as x > 0. 

Emden’s solutions form a grid as in Case IV. 

Case VI. o+ 2<0, 20+n+3< 0, ot n+1< 0. 

The types and arrangement of the solutions are as in Case V except 


that form (W) is replaced by 
O(x) ~ Dz, (X) 


where D is one of the constants of integration. 
12. Tabulations of the behaviour of the proper solutions of 
the equation 
0” +276" = 0 (n> 1, n= p/q, one of p, g even). (2) 
When 7 is not odd the only possible change in the arguments con- 
cerns the behaviour of a solution which crosses or attempts to cross 
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the axis 2=0. When n= p/q, p even, the solutions can cross the 
axis, but on putting 6 = —¢ (¢ > 0) the solutions on the other side 
of the axis satisfy the equation 

$"—x"$" = 0, (8) 
and can be best considered under that heading. When n= p/q, 
q even, proper solutions cannot cross the axis § = 0. We can therefore 
agree to study the proper solutions of equation (2) with the extra 
restriction that they ultimately remain positive. It has of course 
already been proved in I that no proper solution of (2) for these 
values of » or of (8) for any value of » can oscillate continually as 
x —> 0, so that proper solutions 6 = 6(x) must be ultimately mono- 
tonic and also possess a monotonic 6’(x). The results naturally are 
that positive proper solutions of (2) in this case must have one of 
the forms available for monotonic solutions for p, g odd. It is perhaps 
worth while to repeat the relevant parts of Table I. 


Table II 


The form as x > «© of positive proper solutions of 
6” +276" = 0 (n>1, n= p/q, p or q even). 
Case I. o+ 2>0, 20++3 > 0, o+n+1>0. 


There are no proper solutions remaining positive as x > 00. All solu- 


tions again cross or attempt to cross 6 = 0. 

Case II. o+- 2<0, 20+n+3>0, o+n+1>0. 

There are positive proper solutions of the forms (Q) (special) and (R) 
(Emden’s), but all others again cross or attempt to cross 0 = 0. 

Case III. There are positive proper solutions of the forms (Q) 
(special), (S) (Hmden’s), and (T) (with two arbitrary constants). 

Emden’s solutions form a grid separating the positive monotonic solu- 
tions from the others which must try to have another zero. 

Cases IV-VI. The forms of the solutions are in each case identical 
with the positive forms in Table I. Emden’s solutions form a grid 
separating the positive monotonic solutions from the others which must 
have another zero. 

13. Positive proper solutions of the equation 

6” —x79" = 0 (n> 1). (8) 
It has already been shown in I, by arguments that transfer im- 


mediately to rational n, that the equation (8) cannot have any solu- 
tion oscillating about 6=0. We can therefore confine attention to 
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the region @ > 0 and need not distinguish between the various forms 
of n. The solutions are ultimately monotonic and therefore 6(2) > 0, 
lor co. Further, @’(2) is also monotonic, for it cannot vanish more 
than once unless 6”(2), and therefore 6(~) vanishes between each pair 
of roots. Thus also 6’ — 0, 1, or oo. 

The further analysis in I of the monotonic asymptotic forms is 
rather incomplete, and is based on certain general theorems which 
are not immediately applicable for rational n, though the analogous 
theorems most probably continue to be true. The discussion may, 
however, be conducted throughout quite simply by the use of Hardy’s 
theorem* and the usual simple arguments. We shall therefore give 
here a condensed statement of a complete discussion carried out in 
this way. 


Case I. ot 2>0, 20++3>0, o+n+1>0. 


1 ’ 
Here 0=[w(w—1)|""277,, w=— 


v” —bv’+c(v—v") = 0 (6>0,c>0), 

thie independent variable being t= logz. The v(t) of a proper solu- 
tion of (62) can only continue to have zeros as too, if v’(t) has 
a series of zeros and therefore if v = v(t) oscillates about v= 1. But 
this is impossible from the sign of v"(t) near zeros of v’(t). Therefore 
v = v(t) is ultimately monotonic and v(t) > 0, 1, / (JA0,1) or oo. It is 
easily proved in the manner of § 6 that v(t) > 0, 1 (JA 0, 1) are impos- 
sible; v(t)>1 or v(t)=1 is possible. By writing v’ = y equation 
(62) becomes dy 

——bn+c(v—v") = 0, 

nz, —ont+el ) 

which can be studied by Hardy’s theorem as v > 00. It is easily found 
that the only possible form in which v > oo is 


2 

2 (n+1)\#r3 ~ 
10) — | —— — t,—t) »- P 
W)~ [5 (SE) Go (63) 
which of course is not a proper solution. Thus the only positive 
proper solutions are 1 | 

~ [(o+2)(o-+-n+1)]n-a $5 
ed bee 


* The use of this theorem, however, does not extend easily to equations 
in which x? is replaced by f(x), where f(x) ~ x°. 


(64) 
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All other solutions are such that 6(x) vanishes or 6(x) > oo for finite 
x. A closer examination of the solutions for which v(t) ~ 1 shows 
that they are a family depending on one arbitrary constant only. It 
is easy to prove that there is just one member of this family through any 
point and that it separates the solutions for which 0(x) > oo for finite x 
from those for which 6(x) vanishes for finite x. 

Case I (bis). o+-2=0, 20+n+3 > 0, 0+n+1>0. 

Writing x = e the equation becomes at once 

6” —0’—6" = 0. (65) 

It is easily verified that the nature and arrangement of the solutions 


remain the same except that the form (64) is replaced by 


1 


] n—1 
O~ 1 66 
li — l)log zs (66) 
where C is the single constant of integration allowed to the family 
of proper solutions. The form (63) for the solutions that tend to 
infinity as ¢ + ¢, survives in the form 


9 2 


A(t) ~ Schr (67) 


Cases II, IV. c+ 2<0, 20+n7+3 > <0, 0+n+1>0. 


eS o+2 
Here 0=[o(1—am)]"—2v, a= —— = > 0, v> 90, 
n— 
v” —bv'’—c(v-+v") = 0 (c > 0). 
The sign of 6 is unimportant. Again the nature and arrangement of 
the solutions remains the same as in Case I, except that the family 
of proper solutions depending on one arbitrary constant has the 
asymptotic form 
Arye, 
6(x) = A+ = 
(oF 1)(o-+2) 
Case III. o4+2<0, 20+n+3=0, o+n+1>0. 
The nature and arrangement of the solutions remain the same 
except that (68) is replaced by the exact form 


1 
re k(n+1)C2 ]n-1 
9(x) a less A(n a ? (69) 


C being a constant of integration. 


(a ~ 1). (68) 
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Case V. o+ 2<0, 20++3 <0, o+n+1=0. 
Here 0 = 20, v>0, 
v’+v'’—v" = 0, 
but the results are identical with those of Cases IT and IV. 


Case VI. o+2<0, 20++3<0, o+n+1<0. 


SS 
Here 6=[w(w—1)|""27y wow > lv >, 
v” +bv’+c(v—v") = 0 b,c>0, 
and for the first time there is some variety of solutions possible. 
Study of the equation in v shows that v > 1 monotonically and v > 0 
are both possible, but that if v > oo the form is (63) and the solution 
is not a proper one. The alternative v > 0 leads to 


goin +2 
(o+-n+1)(o+n-+2) 
A, B are two constants of integration, but B is only effectively present 
if o+n+2<0. The alternative v > 1 (with v = 1) leads of course to 
1 


~ [(o+2)(o+n+1)]2- -2 
00) = [I coma (71) 


O(x) = Ax+ B+a (a ~ 1). (70) 





which represents a family of solutions depending on one constant of 
integration. The member of this family through any point separates 
the solutions through that point which tend to infinity for finite x 
from the solutions of the family (70), some of which, of course, may 
cross or attempt to cross the axis 6=0 and so leave the region 
under discussion. 

We may summarize these results in Table ITI. 


Table III 
The form as x > © of positive proper solutions of 
6” —«°9" = 0 (n = p/q > 1). 
All the solutions are ultimately monotonic. They cease to be positive 


proper solutions either by crossing or attempting to cross the axis 0 = 0 
or by becoming infinite for finite x, in which case they take the asymptotic 


form 2 2 
_at2 9 n+1 4]n-1 Xo ~n-1 


- 
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Positive proper solutions have the following forms: 


Case I. o+2>0, 20+7+3>0, o+-n+1>0. 


— 1 = Lo 
(a) ~ [ot2otnt+D] 4 -S 
(n—1)* 
with one arbitrary constant; curves separate solutions for which 
(a) > 00 from those for which 6(x) > 0 as x > xp. 


Case I (bis). o+-2=0. As for I, except that 


xn J (Z) 


1 


1 n—-1 
(x) ~ lax ioerea 


Cases II-V. o+-2<0,0+n+1>0. As for I, except that 
. Anytt2 
Oz) = A+— ~ 1). 
=!teriein 4 
When 20+-n+3=0, 
1 
4(n-+ 1)C? n—1 
Ox) = |-— . : 
(x) re = 


Case VI. o+-n+1<0. 
1 
| 9 1 a? — _ot2 
(2) = [Ter | ) a 


n—1 


> 


(n—1)? 
with one arbitrary constant, or 

Anyttnt2 
0(2) = Az+ B4+-—— — a~ 1), Q 
*) tepaHierary ||P ) 
with two constants of integration A, B. Curves of the family (Z) separate 
solutions for which 0(x)—> cas x > xy from members of the family (Q). 











A NOTE ON THE HYPERGEOMETRIC AND 
BESSEL’S EQUATIONS (II) 


By 
J. L. BURCHNALL (Durham) and T. W. CHAUNDY (Ozford) 


[Received 26 May 1931] 


1. Introduction 
In part I of this Note* we called attention to certain elementary 
propositions in the theory of the generalized hypergeometric equation 
I](6—a)y = a" TT] (6+b)y, (1) 
which enable us to transform equations of this type, by substitutions 
of the form y = [1(8—c)z, (2) 
into equations of similar or simpler character in z. 
Our arguments were based on the main theorem I (2), namely: 
If z is a solution of the equation 
(5—a) f(d)z = xg(5)z, 


x. 
then y = [] (6—a—sm)z 
s=0 


is a solution of the equation 
(5—a—rm)f(d)y = x” g(8)y. 

In part I we applied this theorem to a special class of equations 

(1), namely, those of the generalized Bessel’s type 
ITI(—a)y = 2"y, (3) 

employing it mainly to throw light on the classification of such 
equations by the types of function entering into their general solution. 

In this part II we undertake a similar inquiry for the more 
general equations (1), in which, however, some of the simplicity of 
the earlier inquiry will be lacking for the following reason. The sub- 
stitution (2) derives the complete solution of the y-equation from the 
complete solution of the z-equation, unless either the order of the 
z-equation is less than that of the y-equation, or else a solution of 
the z-equation is annihilated by the operator [](8—a—sm). 

With equations of the generalized Bessel’s type (3) neither of these 
possibilities occur. On the one hand the order of the y-equation is 

* Quart. J. of Math. (Oxford), 1 (1930), 186-95. References to results therein 
will be given in the form: I (4). 

3695.2 U 
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the difference of the orders of the 5-operators on the two sides of the 
equation, which is invariant for substitutions of the type (2). The 
z-equation cannot therefore be of lower order than the y-equation, 
and in cases of interest is of the same order, when consequently it is 
also of the generalized Bessel’s form (3). Therefore, on the other hand, 
no solution of the z-equation can be annihilated by the operator 
[[(5—a), for this operator annihilates only finite sums of the form 
> A,2*, whereas the generalized Bessel’s equations have no finite 
solutions of this form. 

Thus, as we have found in part I, the complete solution of the 
z-equation leads at once to the complete solution of the y-equation. 
With equations of the more general type (1) this need not always be 
the case, and we must then fill in the ‘missing’ solution (or solutions) 
from other sources. 


2. Amplification of Theorem I (2) 
We begin by extending the theorem I (2) into a form which takes 
account of these missing solutions. 
We evidently lack a solution y, if and only if some z is a finite 


sum of powers atom — 


For this the necessary and sufficient condition is that g(5) have a 
factor 5—a—sm, where s has one of the above values. If there are 
several such factors, suppose s = q is that with greatest s, and con- 
sider the equation 


(8—a—rm)f(8)y = «(5 —a—qm)g,(5)y. (4) 


| aac | 
The substitution y = [] (6—a—sm)z 


s=q 


now gives for z the equation 
(8—a—qm) f(5)z = «”(8—a—qm)g,(8)z 
or sufficiently, if z, = (6—a—qm)z, 
f(8)z, = v7" 91(8)2- 

In this case, then, as compensation for missing a solution, we have 
been able to lower the order of the equation. 

To replace the missing solution we observe that the full equation 
for z, is actually 


IT (8—a—sm){f(8)z.—2"9,(8)} = 0, 
1 


8=q7 
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or sufficiently, 
[5—a—m(q+ 1) {f(8)z,—2™ 94(3)21} = 9, 
which has the first integral 
S82, = 2" gy(8)z, + Cartas, (5) 


where C is some constant. The substitution connecting y, z is now 


gs TI @—a—sm)zy, (6) 


s=qt 
and z, is annihilated by the operator on the right, only if it be a sum 


of powers 
] gatmath), | gatm(r—), 


If such a z, be a solution of (5), then, whether C be zero or not, the 
last of the above powers must be annihilated by g,(5). Thus g(8) 
must have had a factor 5—a—sm with s >q+1, which is contrary 
to our hypothesis about g. Consequently the complete solution of 
(5) with C 40 gives, through (6), the complete solution of (4). 
Replacing g,, 2, by g, 2 we may enunciate this result as follows: 

(1) Jf q, r, 8 be integers such that q<r and g(5) have no factor 
5—a—sm, where q<8<r, then the complete solution of the equation 

(5—a—rm) f(3)y = x"(5—a—qm)g(5)y 
¢—1 


is given by y = [][ (8—a—sm)z,* 


s=qt1 
where z is the complete solution of the equation 
S(8)z = x™ 9(5)z4+- Cattmat) 


and C is a constant other than zero. 

If f(5), g(8) themselves have each a factor of the type 5—a’—r’m, 
5—a’'—q'm, where q’, r’ are subject to the conditions stated in (1) for 
q, r, then we can once again apply the procedure of (1) to diminish 
the order of the equation, adding a further term independent of z on 
the right. The term Cx*+#™+) already on the right is replaced by 
some X defined by an equation 


Il (5—a’—sm)X = Cxatma@ +1), 


s=q+1 


In general it is enough to write 
X = C’xe+ma+), 


* With the usual convention, of course, that a product of no factors is unity. 
U2 
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But, if the operator on the left contain a factor 5—a—m/(q+1), we 
have to take X of the form 

X = Cz*+ma+) log x. 


We can clearly proceed in this way so long as f(5), g(8) possess 
corresponding factors of the requisite forms.* 
We can replace (1) by the more general theorem: 


(2) If gq, r, 8 be integers such that q<r and g(8) have no factor 
d—a—sm, where q<s<r, and if p be any positive integer, then the 
complete solution of the equation 


9=1 
1 [| (6—a—sm—rm)f(5)y = x”™?(8—a—gqm)g(8)y 


s=0 


r—1 
is given by y= [[ (6—a—tm)z, 
t=q+1 


where z is the complete solution of the equation 


p—1 
1 [ (6—a—sm—qm)f(8)z = x™? g(8)z4+- Cattmar?) 


s=1 

and C is a constant other than zero. 

We suppress the proof, since no new principle arises. 
There is a theorem, complementary to (1): 


(3) If r is a positive integer and z any solution of the equation 


(5—a) f(8)z = x™(8—a—rm)g(8)z, 


then y= II (5—a—sm)z 


s=0 
is a solution of the equation 
SB)y = x g(8)y. 

In this case the order of the z-equation exceeds by unity that of 
the y-equation, and one solution z must be annihilated by the operator 
I[(6—a—sm). We see from the form of the z-equation that it has 
a finite solution led by x* and terminated by x*+"™. This is evidently 
the solution annihilated. 

In the same way that (2) is the generalization of (1), so we have 
the generalization of (3): 

* Tf there are several such factors of f(6) mutually congruent to modulus m, 


it is necessary to proceed in order of descending r to avoid confusion and 
possible vitiation of the argument. 
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(4) If r, p are positive integers and z is any solution of the equation 
p-l1 
I] (6—a—sm)f(d)z = 2” (5—a—rm)g(8)z, 


s=0 


then y= II (8—a—sm)z 
s=0 


is a solution of the equation 
p-1 
IT (5—a—rm—sm)f(5)y = x™? g(8)y. 
3. Equations soluble by elementary functions 
In analogy with I (11) we now consider the equation 


(8—a,)(8—ay)...(8—a,, )y = 2"(5+6,)(5+5g)...(8+6,,)y, 
where the order of the 5-operators on each side is equal to the index 
of x” and the differences a,—a,, b,—b, are all integers incongruent 
to modulus m. We may write this equation more conveniently as 


m—1 m—1 
EEG —a—1— ally = tots (7) 


where N,, WN}. are integers of either sign or zero, and a, b are any real 
numbers. By I (4), this is satisfied by some 
[> II (8—e)z, 
where z is a solution of 
m—1 


Ll 8—a—r—mN)z = 2” T] (8-+b+-r-+-mN’)z, (8) 
r=0 


r=0 
and N, N’ are the algebraically least integers of the sets (N,), (N;}). 
By I (6) this equation is satisfied by the solution of 

(8—a—mN)z = 2(6+b+mN')z, 

and therefore (8) has the m independent solutions 

2 = 2+MN(] — wy)-2-0-mMN +N) (w™ = 1). (9) 
Unless a+b+m(N-+N’) be a negative integer or zero, the expressions 
(9) for z are not expressible as finite sums of powers of x and therefore 
are not annihilated by the operator [[(8—c). We can thus assert that 
(5) In general, the equation (7) is soluble by algebraic functions. 

In the exceptional case, when a+b+m(N+N’) is a negative 
integer or zero, the substitution y = xy’ in (7) and a cyclic rearrange- 
ment of factors on the right make it enough to consider, in place 
of (7), the equation 


m—1 


m—1 
TT (6—r—mN,)y = 2 TT B-+7-+-mN yy. (10) 


r=0 














294 J. L. BURCHNALL AND T. W. CHAUNDY 
The exceptional case is now that in which N-+-N’ < 0. If we write 
the right-hand side of (10) in the form 


m—1 


am Tl (5—r+m +m rou 
r=0 


and apply (1), we see that solutions are ‘missing’ only when 

N, > —1—Nj,_,- 
Let us then divide the integers r= 0, 1,...,.m—1 into two groups 
(s), (¢), according as 

N,+N js > 9, NAN},-1 <0. 
By (1) and the subsequent discussion it appears that the complete 
solution of (10) is given by 
y=Il [1 (—s—kmy, 

8s k=—Nin-s 

where z is the complete solution of 


T[(8—t—mN))z = x™ TT (8—t-+m+mNj,_ z+ ¥ Cx-s-mNin-s, (11) 
t 8 


t 
No logarithms appear, since the factors 5—s—mWN, that have been 
removed from the right are mutually incongruent to modulus m. 
Operation on (11) with the operator 


—Nn-1-1 
II I] (8—t—km) 
t k=Ni+1 


gives the equation 


, 





a 
(l—a™)JT [I (8—t—km)z= > Cyx-8-Nin-s, 
t k=, 8 
that is to say, ‘ 
C.z-t-™Nm—s 
I] (8—t—km)z = > ———- (12) 


No terms have been lost on the right, since s¢¢(modm). The 
complete solution of (12) thus includes the complete solution of (11). 
In (12) the ‘complementary function, is a sum of powers of x; a 
‘particular integral’ evidently introduces, besides powers of x, the 
logarithms log(1—wax), where w” = 1. Terms involving log2 could 
also conceivably appear in the particular integral. In this case z must 
include aterm al tkm log x, 

which gives rise to a term 2‘+*™ on the right of (12). But there are 
no such terms, since ¢, m—s are never congruent to modulus m. Thus 
no logarithms appear except those of (l—wza), and the complete 
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solution of (12) involves at most these and a finite set of powers of 
x. The same is true of (7) in the exceptional case, and we may thus 
sum up and say 
(6) The differential equation 

m—1 m—1 

I] (6—a—r—mN,)y = 2” TJ (6+64+7r+mN}), 

r=0 r=0 
where N,, Nj. are integers of either sign or zero, and a, b are any real 
numbers, is, in general, soluble in algebraic functions, and in any case 
is soluble in algebraic functions together with the logarithms log(1—wz2), 
where w™ = 

4. Legendre’s equation 
In illustration we may apply the preceding theory to a discussion 

of Legendre’s equation 


8(5—1)y = 2*(8—n)(5+n-+ ly (13) 
in the case in which n is a positive integer. If in (2) we write 
a=0=rF, m=l1, p=2, q = —(n+1), 
f6)=1, — g(8) = (6—n), 


we obtain the solutions 


n n n —s n r dt 
D"(a?—1)", D [(2* 1) lear}: 


recognizable* as constant multiples of P,(x), Q,,(2). 

Since n, —(n+1) are of opposite parity, Theorem (6) is applicable 
to Legendre’s equation, which, as is otherwise known, is soluble in 
terms of algebraic functions together with the logarithms log(1+-2). 
Various formulae} expressing Q,,(x) more or less explicitly in alge- 
braic and logarithmic, parts may be obtained by transforming 
Legendre’s equation into the symbolic form 


ea 8(5—1)...(8—n+1) C\x 
Y= *  “§(8—2)...(8—2n) 1—a?” 
and suitably manipulating the fractional operator. 


If we apply (4) to (13), we obtain the known formula 
D"31Q, (x) = 2"nl(1—a2)m+®, 


(14) 





* For Q,(x) compare Whittaker and Watson, Modern Analysis (1920), 
15.31, Ex. 1. The proof there given is essentially the method of the present 
note. 


+ Cf. Whittaker and Watson, loc. cit., p. 333, Ex. 29. 
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determining the numerical factor from a comparison of leading 
coefficients. 


We may obtain certain less usual formulae for P,, Q,, if we 
separate the cases of n even, n odd. Write first n = 2p and in (1) set 
a= —(2p+1), g=08, r=p-+l, 
f(s) =4, g(5) = 8—2p. 


We then have the complete solution of Legendre’s equation, 


p 
y= II (8+ 2p+ 1 — 2s)z ae u(x 1f—))PAv-ly — 7 *P-l(a3D)Pxz, 
s=1 


where z is the complete solution of the equation 


dz = 22(8—2p)z+ CaP +1, 
and so 


io 8) 


z= O'(2®—1)-?4+ O(e®—1)? | x-*»(a®—1)-P1 dex, 


4 


From the coefficient of C’ we obtain, on comparison of leading coefti- 
cients, the formulae 


2?p!P»,,(x) ae x(a-1D)Pa2P-1(y2 -1)?=2 “2p (a3 D)Px(a?— 1)”, 
which may also be written as 
p!P5,,(x4) = x! Dr{xr-4*(x—1)?}, 
p' Psy (a-*) = 2? +t DP {a-P-4(2—1)?}. 
We may express these in the form of contour-integrals as 
xt f 2?-4(2—1)P dz 
x zP-A(z dz 
Py,(x!) = [ ’ 


Qari J (z—a)P tt 
FA 


xPtt ( z—P-h(z—1)P dz 
refer | eae 
277 . (z—ax)P* 
where in both integrals I‘ is a contour in the z-plane which encloses 


z=x but not z=0. These contour-integrals are evidently special 
cases of those which represent the hypergeometric function. 


Again from the coefficient of C in (16) we similarly obtain 
(—)P22?-14 (2p) | 7 

Qo(x#) = ; LDryr-\(~—1)P | t-P-4(t—1)-?4 

Qap( at) = EAP whoa Wa—1p | £--H(0—1)-9- a, 


x 
92 


a PN (2p)! 4 +P DPy-t-P(1 — r)P oe dt 
p\(4p)! 


0 
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If nm still equals 2p and we apply (3) to (13), we see that, if y is 
any solution of Legendre’s equation of order 2p, then 


z= 11 6—2ayy (17) 


is a solution of (8—1)z = x9(8+4 2p+ 1)z. (18) 


Since (18) is satisfied by , _ x(a2—1)-P-1 
and P. 


2, is annihilated by the operator on the right of (17), we have, 
on comparison of leading coefficients, 


11 (8—28)Qap(2) = (—P2”pla(e?—1)-?>, 


which may be written in the equivalent forms 
D?+1Q4 (2!) = (—)P14ple-PA(a—1)-P, 
DP? Qy,, (a4) = (—)P*14plaP4(e—1)-P, 
Corresponding results for n = 2p+-1 are 
PP oy 43(2*) = D?{ax? 4(a—1)}, 
P! Pop s3(a*) = 2? DP fa? (x— 1)?}, 
Qepa(xt) = TE Deer i(e— 1p | x-P-H(x—1)-P-1 da, 
KP + WF yp 1 Dpy-P-*(1—a)P | x*p+t(1|—x)-P— da, 
p\(4p+2)! : 
DP +2-4Q2y 44(0#) = Epla-P-(1—2z)-P 4, 
DP +P +4 Qa, 4(a-+) — dplar+4(1—ax)--1, 
In similar fashion we can obtain formulae for the associated 
Legendre functions by considering the differential equation 
5(8—1)y = 2°(5-+-m—n)(5-+-m-+n-+ l)y 
satisfied by (1—a?)-*™@P™(z). 


Qepsi(e) = 








SOME SELF-RECIPROCAL FUNCTIONS 
By G. N. WATSON (Birmingham) 
[Received 6 July 1931] 


1. Tuts note contains a few minor developments of a recent paper* 
entitled ‘Self-reciprocal functions’ by Hardy and Titchmarsh. The 
problems which I discuss have arisen out of a question which was 
put to me by Mr. E. G. Phillips; he had encountered the formula+ 


425 J, (nx) , (1) 
n=1 pe pecs 8 ar? 
where p is the greatest integer such that 27p<2 and the dash 
indicates that the first term of the sum is to be halved; and his 
question was whether an analogous formula existed for the series 
obtained from (1) by changing J,(nx) into K,(nz). 
In this paper I first consider the more general series 


2) 


> (4nz)’K,(nz) 


n=1 
and I transform it into an infinite series which is a natural modifica- 
tion of the sum on the right in (1); it is then a trivial matter to 
obtain the transformation in the special case v= 0 by a limiting 
process. It is, however, more interesting to observe that a slight 
modification of the series (2) yields a function F(z) which is self- 
reciprocal in a certain Hankel transform, i.e. such that 


4) 


=! J,(yz)F (y) dy. 


This result suggested to me a further problem, namely, whether, 


if a series f(x) +f(2x)+f(8a)+ 
is self-reciprocal in a Hankel transform, it is possible to construct 
a kernel w(x) for which the individual function f(x) is self-reciprocal, 
so that a 

fix) = | w(ay)fly) dy 


| 


* Quart. J. of Math. (Oxford), 1 (1930), 196-231. 
+ This formula is given by Nielsen, Handbuch der Theorie der Cylinder- ° 
funktionen (1904), 336, with its extension to functions of order v.° 
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I give a formal solution of this problem and of the corresponding 
problem for the series 


f(x) —f(3x)+f(5a)—...; 
and lastly I investigate a transform which I believe to be new. 
2. We proceed to the discussion of the series (2) in which it is 
supposed that z and v are complex with R(z) > 0, R(v) > 0. 
We take Poisson’s summation formula* in a form applicable to an 
even function f(t), namely, 


f(0) +2 > s(n) =2[ fi dt+4 > [ fl0)cos 2ant dt, (3) 


and, for the specified domains of v and z, we may define f(t) for 
positive values of ¢t by the formula 


S(t) = (}tz)’K, (2). 
On the left of (3) we have to interpret f(0) as lim f(t) = $T(v). 
To evaluate the integrals on the right, we emik that 
cos at = v(}zat) J_,(at), 
and that, by a known formula,t 


fala te) yp Ue Pt bt PPG P-W+D 
tr QA+ge-AHD (yy + 1) 


p—A+vr+1 p—A—v+1, é _@ 
<F( on ? 9 — ;p+1; 2 ° 
Hence, by taking A= —v—} and p = —}, we get 


po Pyro +4)™ 
| (4tz)”K, (tz)cos at dt = e+e) 





ny 
0 


0 


We now deduce immediately from (3) that 
JP (v)+2 E (dnz)"K, (nz) 


n=1 = 
1 1 
= P(4)l(v+4)2? Ez +3 » areal ei 


provided that R(z) > 0, R(v) > 0. 
[I remark that I owe this proof to the kindness of a referee. My 


* For a proof of the formula under conditions which are fulfilled in the 
example which we are considering, see J. R. Wilton, Journal London Math. 
Soc. 5 (1930), 276-9. 

+ See my Theory of Bessel Functions (1922), p. 410, formula (1). 
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original proof consisted in taking the formula 


M r 1 . —na 
(inz)’K, (nz) = ra | a 1)+# 


1 
transforming the integral on the right into a loop-integral valid if 
R(v)>0, and so expressing the left-hand side of (4) as a contour 
integral which is evaluated by considering poles on the imaginary 
axis; my own proof was longer than the proof just given. ] 
3. We now obtain a modification of (4) valid when the restriction 
on v is removed. Let M be any positive integer and consider 


tenet" 4)z2” x 


a) M-—1 


es S oxo S(T hase 
xXi= ~ oan 72) y+ m . (2nzr)2”+2m+1f . 


n=1 m=0 * 


> (7) ) 22mC(2y+-2 zm | 


(2a)?”+ 2m+1 





m=0 °* 
Since the general term of the infinite series is O(n-®”-?“-1) when n is 
large, the expression, qua function of v, is analytic throughout the 
domain R(v) > —M with the possible exception of the poles of '(v+- 4) 
and of the zeta functions; and, in the more restricted domain in which 
R(v) > 0, the expression reduces to the expression on the right of (4). 

Consequently, by the theory of analytic continuation, we have 
x 

AP (v)-+2 ¥ (4nz)’K, (2) = PT + he” x 


n=1 
) M-1, 


vy 1 9 > { 1 * — eo gam ) , 
A g2v+1 | ie : \ (22+ 4222)" +4 d A m (2n7)?”+ -2m+1 ig 
M-—1 
LoS (YB) amv +2m+1)] 
ig Zz | m (Qn yer 2m+1 Ss” (5) 
m=0 / 


not only when R(v) > 0 but also when R(v) > —M, except when v 
has any of the values 0, —}$, —1, —3,..., —M-+4; and it is, of course, 
supposed that R(z) >0 

We now consider what happens when v= 0. In (5) take MU=1 
and proceed to the limit. We get 


a 1 of l 1 | 
2> K,(nz)=7i-+2 ; a 
2, So(nz) {i+ 2. \v sees: saa} | ; 


oh ae Peart), 
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and hence 
a 


23K = oy . cose ad . = Ae 
2 > Koln) es E sa | Ferra sea) | - 
+y+log(3z) —log(27), (6) 


where y is Euler’s constant and R(z) > 0. 
It seems unnecessary to work out any of the other limiting cases 
in which v > —m (m= 1, 2, 3.,...). 


4. A comparison of formula (4) with the result* 


(ia d a Oe nfl 


wept TeFh 


(and the reciprocal of this result in the Hankel transform) inevitably 
suggests the construction of the function reciprocal to the function 
defined by the series on either side of (4). Accordingly, let 


F(z) = Ww)+2 ¥ (dne)"K, (nz) 


1 
=1P(d)r +e" | eat? > ar 


where R(z) > 0 and R(v) > 0; and consider 


z) 
J 


V(ay)J,(xy)(ay)"- Bay) dy — («> 0, a> 0). 


We poe to examine the permissibility of substituting the 
second series in (8) for F,(ay) and then integrating term by term. 
We first tak 

— >see —-i< ee oe. 


1 p+ 
Then ye *Ju(xy) a l< lye ey) | 
| (a®y?-+- 4027 2)+h) (2nz) Hers) ? 





and | ¥ Iuley) — ly"? J(xy)| 
é | (a2y?+- 4n2772)" : i) —_ ak+2)( In) RAv—H-) 


1 ao 
Sj _ | +L ° d 1 
ince | yd, (ay)| dy anc 2 (2nz 


0 


1 
R@vii) 8Fe convergent, we have 


) 


* This result and its reciprocal are given in my Theory of Bessel Functions 
(1922), 434 and 410. 
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ao 


V(ay)J,,( (xy) (ay) yt {> (a2 Ti 272)v + i) 


\m 


S (ay)e+4 
=>. fey uD) a2 am 22 )ytk dy, 


n=1 
0 =) 


from the first of the inequalities; also, since ly*J,,(zy) dy| and 
, 1 
] 

2, Qnn)Re-——D 
n=1 
inequalities that the order of the summation and of the integration 
from 1 to co may be interchanged. 

Combining these results, we get 


2) 
e 


are convergent, it follows from the second of the 


PWLO+)) 4 


V(xy)J,,(zy)(ay)h-*”* HRI (ay) — ay jy" 


ae r 21°(3)P'\(v+3)(ay)4* 
os [x (ry) J,(@y) (a2y2+-4n272)rtt “dy 
so that, by (7), 


io a) 
“ 


Ney. (ey)(ayye-24|B (ay) — POP OT)| dy 


ay J 


TC a ae ee Ee -\ p—v+4 6 -\ 
. ey >> mae caer ssa 
a \2a a _— 


m=1 °¢* 

Now the series on the right is an analytic function of both 
p and v in any finite domains of values of » and v; and by dividing 
the range of integration on the left into the ranges (0,1) and (1,00) 
and using the second and first expressions in (8) in the respective 
ranges, we see that each of the integrals is an analytic function of 
both p» and v in the domains specified by the inequalities 

R(2v) > R(u+4) > Riv) > 0. (10) 
Also 


PMO +3) | Vay), (ay)(aye-?” 

and so : 

P 2 / om \ 2v—p+h Dax 

[ eeu), (ey) ay AR (ay) dy = “() Fal 7) = 


0 
throughout the domains specified by (10). 
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The only point of the inequality R(2v) > R(u+}4) in (10) is now 

to secure the convergence of the integral at the upper limit in con- 
sequence of the term }I°(v) in (8). Since 


SPW) | Vey), (ay)(ay)e* dy 


exists as a ‘generalized integral’ written 
yg 2) 
G | =lim |e, 
A++05 


we find that, if R(u+4) > Rv) > 0, (10a) 


then 


( : (2 o-9-4 2 
G | (ay), (ry)(ay-? Fay) ay = “(5 F,-»i(). 


0 
(lla) 
The results (11) and (11a) show that the functions 
(2m) (x \-#-4 2x 
»\p—2v +4 v(2 ae F ed 
(ax) F,(az), a (5) ri a 
are reciprocal functions in the (possibly generalized) Hankel trans- 
form of order p. 
In particular, if we take » = 2v—} and a= (27), we see that 
when R(v) > 0 the function F(xv(27)), that is to say 


40 (v)+2 >) {nx (4n)}K,(nav(2z)), 
n=1 
is self-reciprocal in the generalized Hankel transform of order 2v—}. 


5. We now give a purely formal discussion of the general problem 
concerning the integral equation satisfied by f(x) when the series 


fla) +f(2x)+f(8e)+... 
is a self-reciprocal function in the Hankel transform of order v. 
We write 


c+at 


bs) =[e-Yle)de, fe) = 55 | $(8)x-* ds, 
and then c—xi 


= ie) ins) 
a) “ 


eS fine) de =F f winx) de =F n-* f ale) de = Le G(s). 


= = n=1 
9 n=1 n 15 0 
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The Hardy-Titchmarsh condition* for s f(nx) to be self-reciprocal 


n=1 
is then that 
hai L(s) $(8) 
2D(3s+ + 3v+4) 
should be an even function of s—}. Since I'(4s)z-8Z(s) is an even 


function of s—4, we have 
(jm)b(s) (477)0-4(1—8) 


Pa Gsth+H TE—b bth +p 














T'(4s)'(4s+-4v4+}) 
1,\}-s_ NS = i = 
Hence " 
4 Of pi Pet ie a 
fe=s5 | rg prbe pepe & 
_1 fay Pe Gsthtpds ~ 
=a | GU DG jorcietbtp | emavo ay 
c—ni 0 
= | w(zy)fly) dy, 
where \ oii 
P I'(4s)l'\(4s+4v+}) 
ry) = L)s 2 2 4 1 -8 12 
am(xy) Qn . (377) PG—tsl(—iet +3) mary) 8 ( ) 


We now have to evaluate the integral on the right. If it is taken 


for granted that the two sets of poles 
8 = —2n, 8 = —v—4—2n (na = 0, 1, 2, 3....) 


lie on the left of the contour, we find, by calculating the residues at 
these poles, that 
2!1(4v+4) 
eo nee . 3 3 ae oe ee ae 
a (xy) = i 1 re oF3(4, vu +2, —hr4 4347 y”) tT 
(9v-+-4) 
¢ } 7 1 y 1 
_ (27)! T(—iv—F) (bray) 2 1,.18-3-8-e08 
+ ; F,(v+1, $v+32, v-+4; 4or*2*y?). 
I D(4v j 3)D(v i 1) 0“ 3 2 4 
2° 4 a 3 
By a formula due to Hanumanta Rao} the expression on the right 
may be replaced by an integral containing Bessel functions, so that 


io) 


aay) = 2 | vet, (etjeos( “7) 2. (13) 


0 
* loc. cit., 198 and 208. 
+ Messenger of Math., 47 (1918), 134-7. 
(1922), 437. 


See my Theory of Bessel Functions 
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Now that this result is discovered, it may be derived from (12) more 
quickly by observing that, when w(x) is defined by (12), we have 
| \(at)J,(at)ar(t) dt 


0 
o ¢c+ni 


1 P(be)P(de+- +2) ie 
~s] J I(t) na yeT(—Iet tye oe 


0 c-—axi 


"G) ran *o =(e)- 


and therefore w(x) is the Hankel transform of 2 = cos". Consequently 


w(x) = 2 [> V(at)J, (2t008( = = 
0 
whence (13) is obvious. 


Therefore, if 5 f(nz) is a self-reciprocal function, f(x) satisfies the 
n=1 


integral equation 


fa) =2 a) fs (ety cxteos( =) fy) F dy (14) 


It may be remarked that 
c+at 


1 —s 
ani | Sei#leye-* de 


is not necessarily equal to > f(nz), as has been taken for granted, 
n=1 


because the pole of {(s) at s = 1 may have the effect of adding an 
extra term to the series; and also modifications of the series like 
those introduced in (3) may be essential since it may be necessary 
to alter the value of c in Mellin’s formula in order to make 

oe 1 c+at 

—; n-*(s)a-* ds 

Daa J =H 

convergent. 
6. We next consider the integral equation which is satisfied by 

f(x) when the series 


S(x)—f(8x) +f (5x)—f(7x)+ 


3695.2 x 
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is a self-reciprocal function in the Hankel transform of order v. As 
before, write 
$0) = fof) 
0 
and then the condition to be satisfied by ¢(s) is that 
L(s)$(s) 
QT (s+ v4) 
should be an even function of s—}4, where, as usual, 
1 1 1 


L@)=7-s+5-at~ 


From the functional equation for L(s) it follows that 
I'(4s-+4)(4/7)*L(s) 
is an even function of s—} and therefore 
8\s-§ T(3s+43)P (4s 
#0)= (7) rae 
Hence, as in § 5, jn 


f(x) = | m(zy)f) dy, 


0 


iy ¢(1—s). 


Vv 
1 
3V 


t) 
s+iv+}) 


where 


cet 


1 7g\e+ Dde-+ 0k 
@, (ty) = -— [ (: ) (e+) 3) (srs j 


a) =T(—4s)P(—4s 


To calculate the value of w,(x) observe that 


fseozcominat f (7rd 


0 c— at 


1 rat 
= ~—sin{—}. 
EF 22, 


= P 1. (a 
Therefore w,(x) is the Hankel transform of —sin( > , and so 
x 2a 


wo 


w,(xy) = [ V(xt)J,(xt)s sin( 3) ©. - (15) 


0 
Therefore, if f(x)—f(3x)+f(5x)—... is a self-reciprocal function, f(x) 
formally satisfies the integral equation 


t)= {| vet fotysin( 2) fy) ) as dy. (16) 
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7. An interesting example* is obtained by supposing that 
> 1(n)f(nx) 
n=1 


is self-reciprocal, where r(n) denotes the number of distinct repre- 
sentations of n as the sum of two squares. 


a r(n 
Since +2 (m) _ 


n=1 


we find that f(x) now satisfies the integral equation 


fla) = | welay)fly) dy, 
where : 


w(x) = 





c+at 
1 of (2) _TP@_ Tast+h+h 1g, 
27 | (=) I(1—s) h—hthtD” 


It follows that 


2) 


J (at) J,(2ct)ar9(t) waz | rea" lds 


7 27? 
=a5"(,/F): 


t t 


a) — 
and hence ws(xy) =3 | Met (et ad dt 
0 


8. We now return to the kernels introduced in § 5 and §6. By 
taking Hanumanta Rao’s formula 


2) 


[ te-AJ,(at)J,(b/t) dt 

: abe T\(bv-+-4p— Iu) 

~ 22a-p HAY (u+1) ogi ogi 

p—v—p , , ptv—p, ,,a% 
ou, — +1, 9 +1; v6 )+ 
a”b+Pl'(3u—}p— dv) 

aD WI Guth tie tl) 

x orp ete 4 eee Ls) 





x 





+; 


16 


* Suggested by example 9 (p. 212) of the Hardy-Titchmarsh paper. 
x2 
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and using Barnes’s asymptotic expansions* for the generalized hyper- 
geometric functions on the right, it may be verified by some rather 
tedious analysis that, when a and b are positive and their product 
is large, = 
[ #-J,(at)J,(b/t) dt = O(a~4eb-1H), 
0 
provided that py, v, p are such that 
— R(v+3) < R(p) < R(u+}). 
and so, when z is large and positive, 
w(x) = O(a), w(x) = O(a). 

Of course, each of the terms on the right in Hanumanta Rao’s 

formula is exponentially large like 
a—i-tph—i+tpe2v(ab) 

when ab is large and positive, but the asymptotic series which have 
this positive exponential as a factor happen to cancel in the com- 
bination of the two functions )/', which is under discussion. 

Again it is easy to see that, when x is small and positive, 

w(x) = O(1)+O(a"*4), a(x) = O(n) + O(@"*4). 

It follows that, when R(v+4) > 0, there is no necessity to interpret 

the integral equations of § 5 and § 6, 


oo wo 


fe) = | wlay)fly) dy, — fw) = | way fy) dy 


0 


in any Pickwickian sense, as would have been the case if the two 
series containing positive exponentials had not cancelled. 


9. In conclusion I remark that Hanumanta Rao’s integral with 
p=0, but pw not equal to +4, as has hitherto been the case, can be 
regarded as a kernel which gives rise to a transform, at any rate 
formally. This transform has not, I think, been noticed previously 
and the kernel is a generalized hypergeometric function of a type 
quite different from the kernels discussed by Fox.t 

Let R(u+3)>0 and Rwv+})>0 and define w,,(xy) by the 


integral 3 
7 (Ay) at 
Wy (ey) = A(ay)! | Slat, ( *) 7 
0 
where A is any positive constant. 


* Proc. London Math. Soc. (2) 5 (1907), 59-116. 
+ Ibid. (2) 29 (1929), 401-52. 
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It is obvious that 


w,, (ty) = A(ay)! f J,(0)J (2) 


Let f(x) be an arbitrary Suita and let g(z) be its transform with 
the kernel @w,, (ty), 80 that 


g(x) = je pvlty)f(y) dy. 


Then, assuming that the various changes in the order of the integra- 
tions are permissible, we have 


60 


| mleygly) dy 


0 


dt 
- 


=| i (ry), (yz) f(z) dedy 


DODD 


-[fffmmmmnir a Art 


0000 


Now, if we replace y by y7'/A, we have 


j i (<2), (“<) (“*) J, (4) fle) dedy 
: “3 fewys (ue, (yz) f(z) dedy = (7): 


by Hankel’s Sasi. Hence 


[a By rl(ry)g(y) dy = { i (¢T')+J, corg(=) A" dtdT 


0 00 
ow 


(wed (way (P4( 2) Tam 


FF (uc), (we) f (v) dvdu 


Se rhins Se 
<5 © 


Ro 


— 


again by Hankel’s re 

Consequently f(x) and g(x) are (formally) reciprocal functions in 
the transform with kernel By,,y(XY). 

The formulation of conditions with which the formal analysis is 
justifiable is a matter for further investigation. 

















“ ON A TAUBERIAN THEOREM CONCERNING 
DIRICHLET’S SERIES WITH POSITIVE 
COEFFICIENTS 


By K. ANANDA-RAU (Madras) 
[Received 15 May 1931; in revised form 27 July 1931] 


LET «, a, Q,..-, &,, be real constants not all zero, and let the first of 
them which is different from zero be greater than 0. Let us write 
for shortness 1,u = log u, lu = log(l,w),..., l.u = log(l,_,u),...; and let 
L(u) = (1,u)™(/,u)™...(L,,u)™. From the conditions imposed on the 
«’s it is clear that w*L(u) > co as u—> oo. The following is a known 
theorem : 


THEOREM A. Suppose that a, >0 for n = 1, 2, 3,..., and that 


n= 


f(s) = >) ae ae 


is a Dirichlet’s series convergent for s > 0; suppose further that when 


s—>0O A 1 
8 8, 
where A >0 is a constant. Then as n > o, 
AML (A, 
A, =4a,+a,+...+a,~ ra (2) 


The above theorem, with the restriction ,,,, ~~A,,, is due to Hardy 
and Littlewood.* It has been recently shown that it is not necessary 
to include any explicit restriction on the A’s in the hypothesis of the 
theorem.+ The object of the present note is to deduce Theorem A in 
its general form for any Dirichlet’s series from the particular case of 
the theorem for power series, that is, the case when A,, = n. 


* Messenger of Math., 43 (1914), 134-47 (141-6). In an earlier paper Hardy 
and Littlewood gave a proof of the particular case of Theorem A for power 
series: see Proc. London Math. Soc. (2) 13 (1913), 174-91 (180-5). 

+ O. Szdsz, ‘ Verallgemeinerung und neuer Beweis einiger Satze Tauberscher 
Art’, Sitzwngsb. der Bayer. Akad. (Math. Naturwiss. Abt.), 1929, pp. 325-40. 

J. Karamata, ‘Neuer Beweis und Verallgemeinerung einiger Tauberian- 
Satze’, Math. Zeitschrift, 33 (1931), 294-9. 

K. Ananda-Rau, ‘On Dirichlet’s series with positive coefficients’, Rend. 
di Palermo, 54 (1930), 455-61. 

I may take this opportunity of pointing out that the proof given in my 
paper in the Palermo Rendiconti referred to above can be easily adapted to 
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To do this let p = p(n), ¢ = q(n) be defined by 
p—1<A, <p, (3) 
A, +qS PS Antgty (4) 
so that g>0. Further, let 
¢*“.= 4, (se = 1, 2, 3....), 


n—-1<A,<n 
F(s) = S c,e-*, 
1 
C, = ¢+¢,+...+¢,. 


We have* = g a 
cee < SY ae**<c,en-We, 
n—1<A,<n 


a) «oo o 
> c,e- < SY ae < e Y c,e-", 
n=1 r=1 n=1 


F(s)< f(s) < &F(s), 
from which it follows that when s > 0 
F(s)~f(s), 
and so by (1) F(s)~51(2). 


Now F(s) is a power series in e~* with positive coefficients, and we 
apply the particular case of Theorem A when A, =n. We then get, 
corresponding to (2), 
P g to (2) = Ap*L(p) 
P  T(a+1) 
when p-> oo. Since p— 00 as n> 0, (5) is true when n > oo. 


(5) 


prove the following theorem, which may be considered as completing the 
earlier one. 
THEOREM. Suppose that a, >0 for n = 1, 2, 3,..., and that 


2] 
f(s) = & ayer’ 
1 
is convergent for s >0; suppose further that, as s > 0, 
fs) ~ AL(*), 


where A >0 is a constant and L(u) = (lyu)™ (1,u)%2...(L,,u)%™ tends to infinity 
when u—> 0. Let the first of the numbers a, X%g;...5 Xp which is different from 
zero (and therefore greater than zero) be a. Then, asn—>o, 


UyAnta) ~ U(An)- 
This theorem can also be deduced from Satz 1 of Karamata’s paper referred 


to above. 
* Cf. K. Knopp, ‘Grenzwerte von Dirichletschen Reihen bei der Annéherung 
an die Konvergenzgrenze’: J. fiir Math. 138 (1910), 109-32 (115-16). 
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Now by (4) 

C, = 2% = Anig=Ant(Gniit---+On4q) (6) 
with the understanding that when v is such that ¢ = q(n) = 0, the 
expression (@,,,+...+4,,,) in (6) and in the following is to be 
replaced by zero. Since p—1-—> oo with n, we have, similarly to (5), 

C,_.~ eee e) E=e 
T(a+1) Pa+l) 
and so Cy = C,—C,_, = 0{p*L(p)}. (7) 
On the other hand, in view of (3) and (4), 
.>4,4:+..40 
p-1<A,<p 
and so, by applying (7), 
Onze +4n4, = 0{p*L(p)}. 
Therefore, combining (5), (6), (8), we get 


Y , 
A,, = Ci,— (4, rit. +a, +q) 


Ap*L(p) — AAL(A,,) 
(atl) W(a+l)’ 
on making use of the relation p~A,, obtained from (3). 








ON VAN DER CORPUT’S METHOD AND THE 
ZETA-FUNCTION OF RIEMANN (II) 
By E. C. TITCHMARSH (Ozford) 
[Received 2 October 1931] 
1. Iv my previous paper* (which will be referred to as I) I gave 
a simplified form of the proof of the theorem of van der Corput and 
Koksma that, on the lines o = 1—(q+1)/(4Q—2), where Q = 22-1, 
g=2 275 a Se 
{(o-+it) = O(t"? “log t). (1) 
In particular, taking g = 2, we have 
C(t) = O(¢* log t). (2) 
It is, however, known that 
C(g+it) =O), 0<}, (3) 
—in fact van der Corput and Walfisz} proved (3) with 6= 3%. In 
the present paper we give a simpler proof of (3); and, corresponding 
to van der Corput’s later result for the divisor problem, we obtain 
a slightly smaller value of 0, viz. 6 = 7%. 


2. The problem is to obtain suitable inequalities for sums of the 


form > e2mi fin), 


where f(n) is a real function. There are two processes by which such 
a sum may be made to depend on another sum, of the same general 
form but with a different function instead of f(n). One process, due 
in its original form to Weyl, is Lemma 3 of I. The other is an 
‘approximate functional equation’, and is Theorem 1 of van der 
Corput (2). Let us denote these two processes by A and B re- 
spectively. Lemmas 1 and 2 of I are an incomplete form of B, say 
B’. The proof of (2) given in I can be expressed by 
AB’, 

i.e. we first transform the given sum by Weyl’s process, and then 
use Lemmas 1 and 2. 

The proof of (3) due to van der Corput and Walfisz would be 
denoted in the same way by 

ABAAABABAB. 


* See above, pp. 161-73. 
+ See the references at the end of the previous paper. 
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The proof of (3) given below would be denoted by 
ABAAAAB’, 


the AAAAB’ being given by Theorem 3 of I with k = 4. 

All that we require to know of process A is given in paper I. For 
the sake of completeness we shall now give the proof of B. We con- 
tinue the system of numbering lemmas and theorems from I. 


3. Lemma 4. Let f(x) be a real sw with a continuous and 
steadily decreasing differential coefficient f(x) in the interval (a,b), and 
let f'(b) = a, f'(a)=B. Then 

b 
> e27if(n) — > [ etnttser-va dx + Oflog(B—a+-2)}. (1) 
axn<b a-1<v<B+1 


The proof is similar to that of Lemma 1. 
We may suppose without loss of generality that 0<a< 1; for 
otherwise let r be the integer such that a—1<r< aq, and let 
g(x) = f(x)—ra, a’ = 9'(b) = a—r, p’ = g' (a) = Br. 
Then (1) is b 
> orm — > ( e2zitax)-vx} dx + Oflog(p’—a’+2)}, 
asnsb a’ —1<v< P+ 9 
i.e. the same formula for g(x); and 0 < «’ <1. 
With the notation of Lemma 1, we have 
b b 
2% 


’ ~ 
2at | x(x)er™ OF" (x) dx = — > — | sin Qvra e?7*/@f" (x) dar 
Vv 


v=1 a 
b b 


bs ] +f , 
: | eit Mx)—va\f f (~)—v} dx + 


-> 1 erritix)+va}f’(q) dx + +> 
vy prt 


v=1 re v 
b b 


a 
a (He)—va} dy 4. > ; e2ritNa)—v2)f" (a2) dx 
v<Bt1 § v> 734 +1 
=34+3.+2,+ S. 
say. Since f’(x)/{f’(x)+v} is monotonic, the second mean-value 
theorem gives 
b 
F(e)_ d e27 ila) +v2} 
Jf (@)+v 


steal = Oflog(8+1)}. 
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Similarly 
f'(2) di e27itKix)-vz} 


f'(x)—v 


SF) = Oflont+1). 


Alo ¥= -S"]- > 0(=) = Oflogis+1). 


271 
v<Bt+i v<Bti 
As in the proof of Lemma 1, the left-hand side of (2) is equal to 
b 
» e27i f(n)__ e27if@) dx +4-O(1). 
asn<b . 
Hence ~ 


e27i f(n) 
b b 
oat f eii@dr+ > i) erritf(x)-vz} dx + Oflog(B+1)}+0(1), 


0<v<B+1 


a<n<b 


a 


and, since 0 < «a < 1, the result follows. 


4. THrorem 4. Let f(x) satisfy the conditions of Lemma 4, and let 
f’ (x) and f’ (x) be continuous, and 


ms <|f"(x)|<Am,, — | f"(x)| < Ams, 
fora<u<b. Let n, be such that f'(n,)=v (a<v<f). Then 


Qri f(n) tat e27i{f(ny)—vny} 
> e27t = ¢ Seerere~ 1-00 y4 


axn<b a<v<p |f’(n,)|* 
+ Oflog{2+(b—a)m,}]-+-O{(b—a)mim}}. (1) 


Starting from the result of Lemma 4, we may replace the sum on 


the right of § 3 (1) by 
b 
| etrtteroas dx 
atl<v<f-1 » 
with error O(m;*); for we omit four terms, each of which is O(m ;?), 
by Lemma 2. The new sum is empty unless B—a > 2, so that we 
now suppose that this condition is satisfied. We write 

b mv-5 = nv t 8 b 

[errr de= [ + [ + [ =Atdtds 


a a ny—-5 = ny +8 
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Then b 
1 e27i{f(x)—va} { l | 
Don [ Se 
Sie 79) 


l ) l 
_o/ a. i 
(5/f"(n,+-68)| J (5m) 
and a similar result holds for J,. Also 


8 —§ 
J, = | e27iif(ny-+)—vny—vt} dt — e27i{f(ny)—vny +407’ (ny)4 407'"(nv +90} dt 


—§ F) 
8 
a e277 f(ny) -vny} | em I'mv)f] + O(Bmg)} dt 
-8 
8 
= e2ritsin»)—vny) [ erie f'n») dt + O(54ms), 
-§ 
and, since f”(n,) < 0, 
re) 1 7e*\f"(ny)! 
emi i (ny) at = - a 


e ti l e-tu e-tzi 1 
= GF [ —— du = ——_ + o( ). 
f’(m,)|* {al f"(n,) |}? oF Vu \f’(n,)|* Ooms, 
770"|f" (ny)! 


Altogether we have 
b 


e27i{f(ny)—vny} hari 1 
[ e27iti(x)- vx} dx = vee ut -" o( 
” 4 
\f (n,) m 


+ O(84ms;). 


dm, 


© 
a 


The two O-terms are of the same order if 5 = (mm )-*; and if this 
value of 6 is such that 
a<n,—sd, : (2) 
we obtain 
t 
. a : e27i{f(ny)—vny}—fnt 
e27its(x)—vax} dx nee lle 
en 14 
7. lf (n,) | > 
a 


If (2) is satisfied for all values of v such that a+1<v< B—1, we 
obtain, on summing with respect to rv, 
b 


[ e27i{i(x)—va} dx 
atl<v<B-1 * 


a 


+ O(msz*tm). 


e27itf(ny)—vny} 


Toye 


= enti +O{(B—a)mztm}}. (4) 


ati<v<p-1 











ON RIEMANN’S ZETA-FUNCTION 317 
The limits of the sum on the right may be replaced by (a, 8), with 
error O(m;*); and 
B—a=f'(a)—f'(b) = Of(b—a)m}, 
so that the last term in (4) is O{(b—a)mim}}. The theorem then 
follows from Lemma 4. 
Suppose next that there are values of v for which n,+85>b. The 
argument proceeds as before, except that we now obtain a term 
b—ny 
| eri 1" (ny) dt 
—3 
instead of the corresponding integral over (—5,5). We therefore have 
to consider an additional error 
) 
| ev) dt = O{mz1(b—n,)“}. 
b-—nv 
If ¢4(x) is the inverse function of f’ (x), so that |xh’(x)| = | f’(x)|-+ > 1/mg, 
then = b—n, = (a) —y(v) = (a—v)'(€) > (—a)/mg. 
Hence the sum of these terms is 


0( > *~.)= Oftog(e—a)} = Oflog(2-+(0—a)m}, 


a 
ati<v<B-1 
and the theorem again follows. A similar argument disposes of the 


case n,—5 <a. 
27 


5. THEOREM. C($+it) = oli), 
Consider the sum 


b b 
7, co > n-t — > e—itlogn (a< b< 2a). 
1 =a n=a 


By Lemma 3 ‘ et | ) 
¥,=0(8)+0((§ >| 221) } 


b-s 
where p<b—a, and } = > ¢ Maginis)—toge), 
n=a 


We now apply Theorem 4 to >. We have 


fle) = —5- flog(e+8)—loga}, f'@)= 5, 


f"(x) = st 3x°-+ 3x88? 


7 x3(a+s8)8 


"ni . _. & 2a-+-8s 
P(e) = Qn x?(a+-s)?’ 
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Since s < p< b—a <a, and x >a, we have 
Asta~* <| f"(x)|< Asta-, | f'"(x)|< Asta-4, 
Theorem 4 therefore gives 
4 ie e27i dv) { g - ( 
<—" D sry any) 7 oMerS 


where ¢(v) = f(n,)—vn,, and 


as 


e st . st 
a J (6—s) — 2a(b—s)b’ B = f (a) _ 
We next consider the sum 
Zz e27i Pv) 
agvey 
The numbers n, are given by 


2na(a+s)" 


st 
= V; 


27n,(n,+8) 


Hence 


; st)! (1 78 1 
¢*(v) = —_—___—____. = (s Bila— a med 


Qnv*(s®-+-2Qst/av)! — 2(27 4t ye 7? 

and if a< Avt, so that s/t << A/vt, it is clear that 

eo. , _ K,(st)? én 

+ <|d*Mr)| < er - (t >t,), 
K,, Kz, and t,, depending on k only. 

We now use Theorem 3. We have 
R, = Of(st)ta-* +} = Ofa*+1(st) 4} 

the constant implied in the O depending on & only; and a similar 


result holds for 1/r;. Hence 
1 


1 ' eS = 
avis) — 9 [St (Ck? \R3] , ol (sty Ratt) k= 
2. ols iF) (+(e) (an) 


yi kak 9) ( a-ty kh Bg Skt) 
= O | (st) 4K—2 q@4K—2 +4 O | (st) K '4K-2 qK iis © 
Also |f’(n,)|-? is monotonic and of the form Of{(st)-*a!}. Hence, by 
partial summation, 

oni div 1 k_ 2k 

e27i $v) 


f'(n, ) fee O{(st)' 4K—2q,4K-2 
B vii 


1 


+Of(st)? * ‘sK-3 gk 
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We now insert this in (2). Observing that the logarithmic term in 
(2) is always of smaller order than one of the other terms, we obtain 
ks 2k+1 k 2 1 ett 

>.= —_ 0 | (st)? 4K— 4K—2 q1K— 2 “4) +0 (st)? rE +ik- 41K—-2 q@K— 2 4K—2 + 
+ Of(st)ta"} + Of(st)ta-t}. (4) 
Finally we insert this in (1). Summing with respect to p and dividing 
by p has the effect of replacing s by p in each term. We then multiply 
by a; and we can take the square root of each term separately, since 

(X+Y+...)§< X#4+ Yi... 


The result is 
2k+1 


ee 
Y, = Opa) +0 (p t) 4 &K—4 q@8K- saa) 


& 1,1 2 


+O| (pt) 2K *8K=a g@K* 4 8K=4) + Of(pt)-tak} + Of(pt)tad}. (5) 


The first two terms on the right of (5) are of the same order if 
6K—2k—-4 — 2K—k-1 
pine [a 6K—k-3 } eK-K-3| (1<p<b—a), 


and they are then of the form 


, _8K-1_-12K—-k-1 
6K—k-—3 426K—k-—3 
a t , 


We have to take a as large as Avét, and they are then of the form 
, 5K—k-2 


Olt 26K —k— 9), 
For k = 1, 2, 3, 4, 5,... the index has the values 
3 5 17 73 
' a. ae 
and of these }7 is the least. We therefore take k = 4. Then 


36 —_11 , il 
on [awe “| (a>t"), 
23 11 147 1338 169 15, 
>= O(a"t*)+0(a 328 i) + O(a 184 ")+0(a at), (6) 
This also holds if p > b—a, since then 
36 _11 
Y, = O(b—a) = O(p) = O(at 4" 
which is of smaller order than the third term on the right of (6). 
It follows from a by si summation = for b < 2a, a > #16, 


b 





I 
; 3 


and we obtain 


a “am (7) 


n=a 
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Applying this with a= N, b=2N—1; a=2N, b=4N—1, 
last value of b being N’, and adding, we obtain, for N > 11/36, 
© tl). o(w-%, 2 #:,-B). o(y-ayt 
77 82 482 y 328, 328) | y 2\ | 4 
pa = O(N’) + 0(N 2) + O(N E 4 o(n *e*). 
n n 
2 (8) 
If N’ < Avét, the first term on the right is Oli). The second 
term is also of this form if N >¢‘/”. The remaining terms are then 
negligible, and we obtain 
1 , 27 
164 
> a= ol). (9) 
ti<n<Ath 
We require a subsidiary argument for the corresponding sum with 
n <t/%, and Theorem 3 is sufficient for this purpose. In fact (1) of 
I, § 4.1, with k = 3, gives 
N’ 
] , 3 a , 4 
714,14 Ty — 20 
Da O(N"") — (w<e'), 
n=N 
> 1 19 i 
“nce — 119 . . 
Hence 2, ht ol: ) (N <Alt ), 


and, summing in the usual way, 


Le ~T o(, Jogi) — - o(t*). 


Hence 


The result now follows from the approximate functional equation* 
of Hardy and Littlewood 


(3+ %t) * 7 aX sz nt at Ot-*); 


N<N (she n<n (sh) 


where x = O(1). 


* See E. C. Titchmarsh, The Zeta-function of Riemann (Cambridge, 1930), 
32-4, for a proof, and references to other proofs. The proof there given 
should be corrected by taking the integral referred to at the top of p. 34 along 
(x, «+et'% 0) instead of (x, x+ic). 
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This Treatise is in the main concerned with the forms and analytical 
properties of the functions which arise in connection with those solu- 
tions of Laplace’s equation which are adapted to the case of particular 
a | boundary problems. The investigations take account of the functions 
which are not, as was the case when they were originally introduced, 
confined to the cases in which the degree and order are integral. It is 
hoped that, although of a purely mathematical complexion, it may be 
found to be of use to Mathematical Physicists who are primarily con- 
cerned with applications. 


| CARTESIAN TENSORS 
By HAROLD JEFFREYS = 
i] Demy 8vo. 55. net 


The structural simplicity of equations in mathematical physics when 

written out in full Cartesian form is often hidden by the labour of 

writing out every term explicitly. Attempts have been made to reduce 
this by vector algebra. But the use of tensor notation with the summa- 

tion convention is as simple, and it is the object of this work to illus- 
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| 

| 

| 

| 

| 














HIGHER COURSE GEOMETRY 
By H. G. FORDER 


In Two Parts. Crown 8vo. Complete 6s. 
Part IV, 25.6d. Part V, 4s. 


The first three parts of this work, comprising A School Geometry, were 
aT published in September 1930, and are sufficient for pupils taking the 
aT Matriculation or an equivalent examination. Parts IV and V, now 
|| issued, are designed chiefly for Sixth forms. A thorough treatment 

is given of Similar Figures, Solid Geometry, Geometrical Conics, and 
Hl the Theory of Inversion, thus obviating the need for separate books 
| on these subjects. Unusual care has been taken with the diagrams. 








|| CAMBRIDGE UNIVERSITY PRESS 












































Printed in Great Britain at the University Press, Oxford, by John Johnson, Printer to the University 


1, 








